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USERS  GUIDE:  SIMULATION  MODEL  FOR 
AMMUNITION  PLANTS;  PREDICTION  OF 
WASTEWATER  CHARACTERISTICS  AND 
IMPACT  OF  REUSE/RECYCLE 


1  INTRODUCTION 

Background 

Army  Ammunition  Plants  (AAP)  manufacture 
explosives  and  propellants  for  use  by  all  U.S.  military 
services.  The  Army  began  to  eliminate  pollutant 
discharge  from  ammunition  plants  in  the  early  1970s. 
As  a  result,  most  AAPs  now  meet  National  Pollutant 
Discharge  Elimination  System  (NPDES)  discharge 
requirements;  however,  production  is  not  even  close 
to  design  capacity  at  any  of  the  plants.  Thus,  it  is  not 
clear  that  wastewater  treatment  facilities  could  success¬ 
fully  treat  a  five-  to  ten-fold  increase  in  pollutant 
loadings  caused  by  capacity  production  to  the  levels 
required  by  NPDES  permits.  Therefore,  the  Army 
must  identify  the  effects  of  increased  levels  of  am¬ 
munition  production  on  the  ultimate  amount  of 
toxic  pollutants  discharged  from  AAP  wastewater 
treatment  facilities. 

Objective 

The  objective  of  this  study  was  to  develop  a  com¬ 
puter  modeling  tool  to  help  Headquarters,  Army 
Materiel  Development  and  Readiness  Command 
(DARCOM)  and  AAP  personnel  predict:  (1)  the 
effluent  quality  of  wastewater  streams  as  a  function  of 
ammunition  production  level,  and  (2)  the  effect  of 
various  water  reuse/recycle  options  of  effluent  quality. 
This  report  documents  the  model  formulation,  soft¬ 
ware  development,  and  user  instructions. 

Approach 

Needed  model  characteristics  were  defined,  and  a 
model  was  designed  to  satisfy  these  requirements 
(Chapter  2).  Software  and  supporting  documentation 
were  then  developed  (Chapter  3). 


2  MODEL  FORMULATION 

Defining  Modal  Characteristics 

The  first  step  in  designing  a  model  that  would  meet 
Army  requirements  was  defining  its  necessary  qualities. 
Research  indicated  that  a  model  combining  the  follow¬ 
ing  characteristics  would  best  meet  these  needs: 


1 .  Ability  to  model  any  type  of  water  use  that  may 
occur  in  AAPs,  including  water  for  nine  manufacturing 
processes,  washdown,  cooling,  and  transporting  fin¬ 
ished  ammunition  through  pipes. 

2.  Ability  to  specify  a  water-using  process  as  many- 
distinct  steps  or  as  a  “black  box.” 

3.  Ability  to  change  flow  interconnections  and  re¬ 
run  the  model  to  assess  the  impact  of  recycle/reuse. 

4.  Creation  of  a  user-friendly  interactive  environ¬ 
ment  to  prompt  the  users  for  necessary  data  inputs 
and  changes. 

5.  Use  of  the  C  programming  language  and  UNIX 
operating  system  for  all  software  developed  in  order  to 
simplify  eventual  transport  to  DARCOM  computers. 

A  model  flexible  enough  to  apply  to  the  variety  of 
AAP  water-using  activities  allows  the  user  to  become 
familiar  with  one  tool  which  can  solve  many  problems. 
Designing  such  a  “general-purpose”  model  means  that 
very  few  assumptions  are  made  about  the  character¬ 
istics  of  any  particular  process.  Instead,  the  user  must 
describe  each  process  train  fully. 

General  Description 

The  model  designed  to  incorporate  these  defined 
characteristics  includes  four  basic  steps: 

1 .  The  user  describes  a  process  train. 

2.  The  user  describes  water  use.  pollutant  gener¬ 
ation,  and  treatment  efficiencies  for  each  process  in 
the  train  as  a  function  of  production  level. 

3.  The  program  calculates  water  flows  or  pollutant 
loadings  throughout  the  system,  including  the  dis¬ 
charge;  this  is  displayed  immediately  and  saved  in  a 
file  called  RESULTS. 

4.  The  user  may  change  the  system  characteristics 
such  as  flow  connections,  treatment  levels,  etc.,  and 
re-run  the  program ;  this  allows  the  analysis  of  various 
water  reuse/recycling  schemes. 

The  program  describes  the  flow  network  with  a  set 
of  simultaneous  linear  equations  whose  solution  is  the 
flow  rates  in  each  pipe.  Each  plant  process,  whether 
manufacturing  or  treatment,  is  referred  to  as  a  node. 
Nodes  can  also  be  pipe  junctions  or  other  situations 
where  water  or  pollutant  flow  rates  are  changed.  The 
flows  between  connected  nodes  are  the  variables  the 
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program  solves  for.  Figures  1  and  2  are  simplified 
flowcharts  that  describe  how  the  plant  model  handles 
water  flows  and  pollutant  loadings. 

The  program  describes  the  system  of  flows  by  means 
of  the  connection  matrix  -  the  variable  “array.”  The 
connection  matrix  is  a  two-dimensional  array,  with 
each  node  represented  by  a  row  and  a  column.  The 
columns  represent  nodes  which  flow  is  "from,”  and 
the  rows  represent  nodes  which  flow  is  “to”;  a  “1”  in 
the  matrix  indicates  connection  between  nodes.  For 
instance,  flow  from  node  3  to  node  5  is  represented  by 
a  “1”  in  column  3,  row  5.  The  connection  matrix 
allows  the  user  to  change  flow  pathways  easily  by 
simply  inserting  or  deleting  "l”s  in  the  matrix.  The 
routine  called  Matrxin  creates  the  connection  matrix. 

Once  the  connection  matrix  is  made,  it  is  used 
by  the  routine  Bildsys  to  generate  the  equations 
describing  the  system.  These  equations  are  represented 
by  their  determinant  (variable  “determ”) and  the  vector 
of  right-side  values  (“rhs”).  (In  matrix  notation,  the 
system  of  equations  is  defined  by  “determn[]  (]  * 
var[]  =  rha[]  ”.) 

Several  kinds  of  equations  describe  the  inflows  and 
outflows  at  each  node,  in  a  manner  analogous  to  the 


use  of  Kirchoffs  laws  to  solve  electrical  networks.  The 
first  set  of  equations,  which  is  extracted  directly  from 
the  connection  matrix,  simply  states  that  the  sum  of 
inflows  to  a  node  equals  the  sum  of  outflows.  When 
there  is  a  node  with  more  than  one  outflow,  such  as 
foi  a  recycle,  more  equations  are  added  to  describe  the 
fraction  of  the  outflow  which  goes  to  each  downstream 
node.  At  nodes  where  there  are  water  inflows  to  the 
system,  an  equation  is  added  which  states  thai  the 
inflows  to  the  node,  including  inflows  from  the  water 
source,  equal  the  water  demanded  by  the  node.  At 
this  point,  a  variable  representing  the  source  How  at 
the  node  is  added. 

After  all  these  equations  are  generated,  there  is 
one  equation  for  each  variable  (pipeflow),  and  the 
equations  can  be  solved.  This  is  done  by  the  computer 
with  a  Gaussian  Elimination  routine  in  the  function 
Solvsys. 

Following  a  more  detailed  description  of  each 
routine  of  the  program,  roughly  in  the  order  in  which 
control  normally  flows.  Figure  3  shows  a  simple 
process  train  which  will  be  used  to  illustrate  the  results 
of  each  routine. 


i 
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Figure  1 .  Simplified  flowchart  of  hydraulic  model. 
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Figure  2.  Simplified  flowchart  of  pollutant  model. 
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Figure  3.  Example  process  train. 
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Routine* 

Matrxin 

Matrxin  builds  the  connection  matrix  from  inter¬ 
active  input.  It  was  designed  to  minimize  input  effort 
by  allowing  the  user  to  describe  nodes  with  alpha¬ 
numeric  names  of  20  or  fewer  characters.  The  user  is 
asked  to  put  in  the  name  of  a  node  and  the  name  of 
the  nodes(s)  to  which  flow  goes.  This  is  done  for  all 
nodes;  the  program  builds  the  matrix  from  this  infor¬ 
mation.  The  information  for  each  node  (its  name  and 
the  names  of  the  nodes  flow  goes  to)  is  stored  in  a 
structure  called  “node.”  “node.ndname”  is  a  character 
string  with  the  node’s  name,  and  "node.ndto[]  ”  is  an 
array  of  character  strings  containing  the  names  of  the 
nodes  flow  goes  to.  The  program  is  presently  dimen¬ 
sioned  for  up  to  29  nodes  and  up  to  five  outflows  for 
each  node.  (The  structure  “node”  is  dimensioned  to 
30,  but  the  0  element  is  not  used  likewise,  the  charac¬ 
ter  arrays  tonode  and  flofrac  are  dimensioned  to  6  but 
only  5  are  used.) 

Flows  to  the  outfall  (or  “sink”)  should  not  be  put 
in  as  nodes.  Matrxin  asks  the  user  which  nodes  go  to 
the  outfall  and  automatically  includes  these  flows  in 
the  flow  network.  Likewise,  flows  into  the  network 
from  the  source,  whether  it  is  a  river  or  water  treat¬ 
ment  facility,  are  automatically  included  for  any  node 
at  which  there  is  a  water  demand  (this  is  done  in 
routine  Bildsys). 

To  avoid  having  to  enter  all  the  node  information 
every  time  the  program  is  run,  an  option  has  been 
added  to  allow  the  user  to  read  all  the  input  values  for 
Matrxin  from  a  file.  Matrxin  returns  a  value  to  the 
main  routine  which  is  the  size  of  the  connection  matrix 
which  is  also  the  number  of  nodes  in  the  system. 
Figure  4  shows  an  example  connection  matrix. 


NODE  1 

NODE  1 

0 

NODE  2 

0 

NODE  3 

1 

NODE  2 

1 

0 

0 

NODE  3 

0 

1 

0 

OUTFALL 

0 

0 

1 

Figure  4.  Connection  matrix  built  by  routine  Matrxin. 


Chgmtx 

The  routine  Chgmtx  allows  the  user  to  modify  the 
connection  matrix  either  to  correct  input  mistakes  or 
to  model  different  flow  routes.  The  routine  asks  the 
user  which  node  to  change;  this  refers  to  changing  the 
outflows.  The  nodes  to  which  flow  is  going  are  printed 
out,  and  the  user  is  asked  which  nodes,  ifany,  to  delete. 
The  user  may  then  add  flow  to  any  existing  node. 

This  routine  may  be  called  immediately  atter  the 
input  routine  or  before  staiting  a  subsequent  run  of  the 
program. 

Bildsys 

Bildsys  builds  the  system  of  simultaneous  equations. 
Because  it  is  also  the  main  function,  it  is  known  in  the 
program  as  “main.”  After  calling  functions  Matrxin 
and  Chgmtx,  Bildsys  uses  the  connection  matrix  to 
write  the  first  group  of  equations. 

The  variables  are  labeled  by  the  program  before  the 
equations  are  generated.  Since  each  “1”  in  the  con¬ 
nection  matrix  represents  a  flow,  there  will  be  one 
variable  for  each  “1.”  To  keep  track  of  them,  the 
program  goes  through  the  matrix  from  row  1,  column 
1,  moving  across  each  row  and  numbering  the  “1  ’s 
encountered  in  ascending  order.  These  numbers  are 
then  referred  to  as  “variable  numbers."  Since  the 
“1”  in  the  matrix  is  replaced  by  its  variable  number, 
the  variable  number  for  a  flow  is  expressed  as 
“matrix [row]  [col] .”  There  will  be  a  row  and  column 
in  the  determinant  for  each  variable  number.  Figure  5 
shows  a  numbered  connection  matrix. 

The  first  group  of  equations  (there  will  be  one  tor 
each  node),  simply  states  that  the  sum  o,  inflows 
equals  the  sum  of  outflows  at  each  node.  This  is  done 


NODE  1 

NODE  1 

0 

NODE  2 

0 

NOOE  3 

( 

NODE  2 

2 

0 

0 

NOOE  3 

0 

3 

0 

OUTFALL 

0 

0 

4 

Figure  5.  Numbered  connection  matrix  (variable 
numbers)  in  routine  Bildsys. 
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by  moving  through  the  connection  matrix  and  at  each 
variable  (element  in  the  matrix  with  a  non-zero  value), 
putting  a  +1  in  the  determinant  at  row  =  row  in 
matrix,  column  =  variable  number,  and  a— 1  at  row  = 
column  in  matrix,  column  =  variable  number.  In  other 
words,  the  program  looks  at  a  flow,  finds  the  node  at 
which  the  flow  is  an  inflow,  and  puts  a  +1  in  the 
equation  for  that  node.  It  then  finds  the  node  for 
which  it  is  an  outflow  and  puts  in  a  —1  in  its  equation. 
For  the  bottom  row  of  the  connection  matrix,  which  is 
for  the  outfall,  only  — l’s  are  added,  since  flow  does 
not  go  to  anywhere  from  the  outfall.  The  result  is  one 
equation  for  each  node.  Within  the  determinant,  the 
row  number  is  equal  to  the  number  of  the  node  the 
equation  is  for;  the  column  numbers  correspond  to 
variable  numbers,  with  one  additional  column  for 
flow  to  the  outfall. 

Bildsys  then  asks  the  user  which  nodes  are  processes 
at  which  there  is  a  water  demand.  These  are  nodes  at 
which  the  process  being  modeled  requires  a  certain 
amount  of  water;  this  amount  is  determined  from  the 
production  level  by  the  function  “findmdl”  (explained 
below).  For  such  nodes,  an  additional  variable  must 
be  added  to  represent  the  amount  of  water  coming 
from  the  source  to  this  node.  An  equation  is  added 
which  states  that  the  sum  of  inflows  to  the  node, 
including  the  source  flow,  equals  the  water  demand 
as  determined  by  Findmdl.  The  model  assumes  that 
if  there  is  a  water  demand  at  a  node  and  if  there  is 
insufficient  inflow  from  upstream  nodes  to  meet  the 
demand,  water  from  the  source  is  added. 

The  equation  is  added  by  making  a  new  row  in  the 
determinant  which  has  a  +1  in  each  column  cor¬ 
responding  to  a  variable  number  in  the  node's  row  of 
the  connection  matrix.  The  inflow  is  then  added  by 
putting  a  +1  in  a  new  column,  which  then  represents 
the  inflow;  the  required  flow  rate  is  added  to  the  right- 
hand  side  vector  (rhs)  for  this  new  row. 

Because  another  inflow  has  been  added  to  the  node, 
each  time  this  procedure  is  used  to  add  source  flow, 
the  program  must  go  back  and  add  a  +1  to  the  node’s 
“inflow  =  outflow”  equation.  This  is  done  by  first 
looking  in  the  node's  column  in  the  connection  matrix 
and  finding  the  first  variable  number  in  the  column. 
In  any  row  of  the  determinant  in  which  there  is  a  —1 
in  the  column  equal  to  that  variable  number,  a  +1  is 
added  to  the  column  corresponding  to  the  inflow. 

The  final  group  of  equations  added  describes  how 
How  is  divided  at  nodes  willi  more  than  one  outflow. 
This  is  necessary  to  provide  enough  equations  to  solve 


the  system,  since  an  extra  outflow  means  an  extra 
variable.  The  equations  needed  at  each  node  include 
one  tor  each  extra  outflow,  or  one  less  than  the 
number  of  outflows.  It  is  assumed  that  a  constant 
fraction  of  the  total  outflow  will  go  each  down¬ 
stream  node. 

The  program  moves  through  the  columns  of  the 
connection  matrix.  Any  column  with  a  non-zero  value 
in  more  than  one  row  has  more  than  one  outflow.  At 
this  point,  the  user  is  asked  for  the  fraction  of  the  out¬ 
flow  going  to  one  of  the  downstream  nodes.  This 
question  is  repeated  so  that  all  the  equations  needed 
can  be  written.  Each  equation  states  that  (1  flofrac  [i)  ) 
times  the  outflow  to  node  i  plus  (-1  *  flofrac [i) ) 
times  all  the  other  outflows  equals  zero  ( flofrac  [i] ) 
equals  the  fraction  of  the  outflow  to  node  i).  This 
equation  divides  the  flow  between  node  i  and  the  other 
outflow  nodes.  Figure  6  is  an  example  of  a  determinant 
built  by  Bildsys. 

Modeling  consumptive  use  of  water  is  not  now 
included  in  the  program,  but  doing  so  would  not  be 
difficult.  Water  consumption  either  as  a  constant  rate, 
as  a  constant  fraction  of  water  use.  or  as  a  function  of 
production  level  can  be  modeled  by  altering  the 
“inflow  =  outflow”  equations.  This  process  is  similar 
to  the  way  pollutant  removal  is  handled,  as  described 
on  p  13. 

At  this  point,  the  set  of  equations  is  complete  and 
Bildsys  calls  function  Solvsys  to  solve  them.  The  above 
description  of  Bildsys  has  assumed  that  water  flows  are 
being  modeled  and  there  are  some  differences  when 
pollutants  are  being  modeled  instead.  Instead  of  asking 
if  a  water  demard  function  exists,  the  program  asks  if 
there  is  pollutant  generation  at  each  node.  If  the 
answer  is  “yes."  function  Findmdl  determines  the 
amount  of  pollutant  added  to  the  system  instead  of  the 
water  demanded.  Instead  of  adding  a  variable  and  an 
equation  describing  water  inflow.  Bildsys  modifies  the 
existing  “inflow  equals  outflow”  equation  to  sa\ 
“inflow  plus  pollutant  generated  equals  outflow  " 
This  is  done  by  changing  the  right  side  of  the  equation 
from  zero  to  the  amount  of  the  pollutant  generated  If 
no  pollutant  is  generated  at  a  node,  the  progiam  asks 
the  user  if  the  node  is  a  waste  treatment  process  where 
the  pollutant  is  removed  If  it  is.  the  program  uses  the 
function  Trtmdl  to  determine  the  percent  temoval  ot 
the  pollutant.  The  "inflow  equals  outflow  "  equation  is 
then  modified  to  "inflow  tunes  removal  equals  out¬ 
flow."  This  is  done  by  multiplying  the  i  I  (inflow) 
values  m  the  equation  by  the  removal  ti.kiiou 
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Figure  6.  Determinant  (system  of  equations)  built  by  routine  Bildsys. 


The  final  section  of  Bildsys  prints  out  the  various 
tlows  after  the  system  of  equations  has  been  solved 
and  then  allows  another  run  of  the  program. 

h'indmdl 

The  routine  Findmdl  returns  the  amount  of  water 
demanded  or  pollutant  generated  at  a  node  when 
called  from  Bildsys.  This  is  done  when  the  program 
searches  a  file  of  existing  process  models,  which  are 
two-dimensional  graphs  of  production  level  (as  percent 
ol  full  production)  versus  water  demanded  or  poilutant 
generated.  The  name  of  each  file  is  the  name  of  the 
constituent. 

The  piogram  first  checks  whether  such  a  file  exists 
by  using  the  UNIX  library  function  “open."  If  “open” 
returns  a  negative  value,  the  file  does  not  exist  and  one 
is  created  with  tile  library  function  "creal."  If  the  file 
docs  exist,  it  is  searched  foi  title  lines  lines  which 
begin  with  the  character  “x."  If  the  (tile  is  the  same  as 
the  node  name  in  question,  the  model  is  read.  It'  no 
model  is  found  with  the  node's  name,  a  new  model 
must  be  put  in.  Thus,  the  user  should  use  the  same 
node  names  to  use  existing  process  models. 


If  an  existing  model  is  found,  it  is  graphed  out  by 
the  function  Graph,  and  the  user  may  replace  the 
model.  If  the  user  does  not  want  to  change  the  existing 
model,  it  interpolates  the  desired  water  demand  or 
pollutant  generated  and  returns  it.  In  all  other  cases 
(a  new  file  created,  no  existing  model  found,  or  ex¬ 
isting  model  to  be  replaced),  a  new  model  is  gener¬ 
ated  and  written  into  the  file,  overwriting  the  existing 
model  if  there  was  one.  The  user  is  asked  io  put  in  the 
production  ievei  and  corresponding  water  demand  or 
pollutant  generated  for  arbitrary  number  of  points, 
up  to  20.  When  a  production  level  of  100  is  entered, 
the  routine  automatically  stops,  calls  the  graphing 
routine,  and  asks  rf  the  model  should  he  changed. 
When  the  user  is  satisfied  with  tile  model,  hitidmui 
interpolates  the  water  demand  of  pollutant  generation 
from  the  production  level  specified  by  the  user. 

An  impoitanl  assumption  built  into  Finuimll  (and 
Titmdl)  is  that  when  a  production  level  is  modeled 
which  is  less  than  the  first  point  in  the  model,  the 
water  demand  is  the  same  as  in  the  first  point  ol  the 
model.  In  other  words,  the  program  does  not  assume 
that  water  demand  goes  to  lexo  as  production  goes  to 


12 


zero.  This  assumption  could  be  changed  with  a  minor 
program  modification  in  the  interpolation  routine. 

Trtmdl 

Trtmdl  is  very  similar  to  Findmdl;  the  difference  is 
that  Trtmdl  is  used  to  handle  modeling  of  wastewater 
treatment  processes,  while  Findmdl  is  used  to  handle 
pollutant  generation.  Instead  of  modeling  water 
demand  or  pollutant  generated  vs.  production  level, 
it  models  removal  efficiency  vs.  flow  rate,  a  function 
which  is  more  traditionally  used  in  wastewater 
engineering  and  modeling.  Once  the  proper  model  is 
found  or  put  in,  exactly  as  in  Findmdl,  a  value  of 
percent  removal  is  returned.  To  do  this,  a  value  of  the 
flow  rate  must  be  obtained.  This  value  is  stored  in 
the  structure  “nodes”  under  the  variable  “watrin.” 
where  its  value  was  assigned  the  last  time  the  program 
was  run  for  water  flows.  This  means  it  is  important  to 
model  water  flows  before  modeling  pollutant  flows 
and  removal. 

Graph 

The  function  Graph  is  a  simple  graphing  routine 
which  uses  two  one-dimensional  arrays  as  input.  It 
generates  50  points  between  zero  and  the  maximum 
value  of  the  x=axis,  then  plots  them  out  on  20  rows 
of  50  columns  each. 

Solvsys 

Solvsys  solves  the  system  of  linear  equations  gen¬ 
erated  by  Bildsys,  using  Gaussian  Elimination  on  the 
determinant  “determ”  and  the  vector  “rhs”  of  right- 
side  constants.  Gaussian  Elimination  requires  that  the 
elements  of  the  main  diagonal  of  the  determinant  be 
non-zero  to  avoid  division  by  zero. 

Solvsys  uses  two  separate  methods  to  insure  non¬ 
zero  values  in  the  diagonal.  First  it  moves  down  the 
rows  and  if  it  finds  a  row  with  a  zero  on  the  diagonal, 
it  looks  below  that  row  for  another  row  with  a  non¬ 
zero  value  in  the  desired  column.  If  it  finds  such  a  row, 
the  rows  are  interchanged.  After  going  through  the 
entire  determinant  in  this  fashion,  most  (but  not 
necessarily  all)  of  the  diagonal  elements  are  non-zero 
and  the  Gaussian  Elimination  routine  starts. 

The  routine  starts  at  row  1,  column  1,  and  looks 
down  the  rows  until  it  finds  a  non-zero  value  in  column 
1 .  When  this  happens,  row  1  is  multiplied  by  the 
appropriate  value  and  added  to  this  row  to  make  the 
value  go  to  zero.  The  program  proceeds  in  this  fashion 
down  through  all  the  columns  for  each  row.  If  a  non¬ 
zero  element  is  found  in  a  column  whose  diagonal 
element  is  still  zero,  the  element’s  row  is  added  to  the 


diagonal’s  row  to  get  a  non-zero  diagonal.  When  this 
process  is  finished,  all  values  below  the  diagonal  of 
the  determinant  are  zero;  Solvsys  then  goes  back, 
starting  with  the  last  variable,  and  evaluates  ail  the 
variables.  A  variable  is  evaluated  by  using  the  row  with 
it  on  the  diagonal  (the  row  number  equal  to  the 
variable’s  column  number);  its  value  is  the  right  side 
of  the  row.  minus  the  sum  of  all  other  variables  to  the 
right  of  the  variable  times  their  coetficients  in  the 
determinant;  this  number  is  divided  bv  the  coefficient 
of  the  variable  in  question.  The  values  of  the  variables 
are  stored  in  the  array  "var[] ."  At  this  point,  control 
returns  to  the  main  routine  Bildsys.  which  prints 
out  the  results  and  writes  them  into  a  file  called 
“RESULTS.” 

Help 

Help  is  called  whenever  the  user  replies  to  a  prompt 
with  the  word  "help.”  Help  includes  some  general 
instructions  on  using  the  program  and  specific  in¬ 
structions  for  each  prompt.  The  instructions  for  each 
prompt  are  based  on  the  information  in  the  user's 
guide  (Chapter  3). 

Model  Verification  and  Field  Application 

Enough  data  has  not  yet  been  obtained  from  any 
single  ammunition  manufacturing  process  to  do  a 
verification  study  of  the  model.  There  are  many 
publications  from  AAP  water  pollution  studies,  but 
there  is  no  complete  set  of  information  on  process 
descriptions,  water  use,  pollutant  generation,  waste- 
water  treatment  capacities,  and  effluent  characteristics, 
all  measured  under  the  same  production  conditions; 
such  data  is  needed  to  make  a  realistic  model  run. 
The  problem  is  compounded  because  much  waste- 
water  treatment  modernization,  including  water  reuse 
and  recycling,  has  been  done  since  plants  stopped 
large-scale  production;  in  fact,  the  present  lack  of 
activity  at  AAPs  will  hamper  any  modeling  effort. 

Conversations  with  Government  engineers  at  AAPs 
have  determined  that  the  input  data  needed  to  run 
this  model  is  available,  though  it  may  have  to  be 
extracted  from  operating  logs.  Unfortunately,  these 
logs  are  not  available  to  non-DARCOM  personnel, 
which  indicates  that  the  program  should  be  used  by 
plant  operators  or  overseers,  since  they  are  most 
familiar  with  plant  operations  and  production  levels. 

One  problem  which  will  be  encountered  when  the 
model  is  used  in  the  field  is  that  some  production 
lines  do  not  really  have  gradations  of  production 
level;  i.e.,  they  are  either  on  or  off.  The  model  can 
handle  this  situation  in  several  ways.  If,  for  instance,  a 
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production  line  is  either  on  or  off  during  each  of  three 
daily  shifts,  the  water  use/pollutant  generation  graphs 
can  be  drawn  as  three-stage  step  functions:  for  up  to 
33  percent  production,  the  pollution  generated  by  8 
hours  (one  shift)  of  continuous  operation  is  used;  for 
up  to  67  percent  production,  the  pollution  generated 
by  16  hours  (two  shifts)  is  used;  and  for  100  percent 
production,  the  pollution  generated  by  24  hours  (three 
shifts)  is  used,  (n  such  cases,  water  flows  will  be 
averaged  over  the  whole  day,  which  is  realistic  if  there 
is  much  detention  time  in  the  wastewater  treatment 
system. 

The  program  can  also  be  used  with  flat  water  use/ 
pollutant  generation  curves,  meaning  that  pollutant 
generation  is  independent  of  production  level,  or  that 
production  level  is  irrelevant.  In  these  cases,  where  the 
production  line  is  either  on  or  off,  the  effluent  is 
readily  predicted  by  using  historical  measurements  of 
pollutant  flow  as  the  (flat)  pollutant  generation  graph. 
When  the  line  is  on,  these  measured  flows  will  be 
routed  into  the  treatment  models  to  predict  the 
effluent;  when  the  line  is  off,  no  effluent  will  be 
predicted. 

Another  problem  in  field  use  will  be  that  in  some 
cases,  data  from  only  one  production  level  will  be 
available.  In  this  case,  the  user  will  have  to  estimate  the 
shape  of  the  curves;  usually  this  will  mean  extrapo¬ 
lating  toward  higher  production  levels.  An  experienced 
ammunition  plant  operator  should  know  to  some 
extent  how  water  use  and  pollution  generation  will 
change  as  production  increases;  for  instance,  the 
operator  may  feel  that  the  amount  of  washwater  will 
remain  constant  if  washdown  operations  are  not 
increased  with  production,  while  the  discharge  of  some 
pollutant  from  the  actual  production  line  may  be 
estimated  to  be  linearly  related  to  production  level. 
This  requirement  for  user  judgment  is  another  reason 
why  plant  operators  should  use  the  program. 

Possible  Model  Modifications 

The  plant  model  could  be  further  developed  to 
include: 

1 .  Built-in  models  of  the  manufacture  of  explosives 
and  propellants  other  than  TNT.  Currently,  such 
models  must  be  input  by  the  user. 

2.  A  method  to  account  for  the  consumptive  use  of 
water  and  reaction  kinetics  of  pollutants.  A  kinetic 
model  would  involve  a  major  effort;  however,  a  model 
for  consumptive  use  of  water  could  be  added  easily. 


3.  A  checking  routine  for  process  trains  involving 
reuse/recycle  to  compare  the  recycled  water  quality 
against  the  water  quality  needed  for  process  influent. 

4.  Environmental  fate  of  toxic  pollutants  m  the 
receiving  stream.  A  number  of  adequate  fate  models 
exist,  but  they  are  not  written  in  the  C  programming 
language. 

5.  Costing  algorithms  and  expected  treatment 
efficiencies  for  various  treatment  technologies  used 
to  treat  ammunition  plant  process  wastes. 

3  USERS  GUIDE 

The  plant  model  is  available  as  Experimental  Profile 
Number  30  of  the  Environmental  Technical  Infor¬ 
mation  System  (ET1S).1  ETIS  can  be  accessed  over 
commercial  telephone  lines  by  almost  any  kind  of 
computer  terminal.  ETIS  is  an  interactive,  answer- 
driven  umbrella  system  which  can  be  used  by  people 
unfamiliar  with  computers. 

Accusing  ETIS 

If  you  decide  to  access  the  plant  model  after  inter¬ 
actively  entering  ETIS,  you  need  only  transfer  into  the 
program.  The  following  paragraphs  provide  instructions 
for  accessing  ETIS  information  by  remote  terminal. 

Once  you  have  acquired  a  log-in  and  a  password 
from  CERL's  Environmental  Division,  you  can  access 
ETIS  by  remote  terminal  by  following  the  directions 
in  CERL  Technical  Reports  N-562  and  N-23  (DA 
Pamphlet  200-2*),  and  CERL  Technical  Report 


1  R.  D.  Webster,  et  at.,  Development  of  the  Environmental 
Technical  Information  System,  Interim  Report  E-52/ADA 
009668  (U.S.  Army  Construction  Engineering  Research 
Laboratory  [CERLI ,  1975). 

2J.  van  Weringh,  J.  Patzer,  R.  Welsh,  and  R.  Webster. 
Computer-Aided  Environmental  Legislative  Data  System 
(CELDS/  User  Manual,  Technical  Report  N-56/ADA061126 
(CERL,  1978). 

3R.  D.  Webster,  R.  A.  Mitchell,  R.  L.  Welsh,  E.  Shannon, 
and  M.  L.  Anderson,  The  Economic  Impact  Forecast  System 
Description  and  User  Instructions.  Technical  Report  N-2 / 
ADA027139  (CERL.  1976). 

4  77te  Economic  Impact  Forecast  System- Description 
and  User  Instructions,  DA  PAM  200-2  (Department  of  the 
Army,  December  1976). 
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N43.5  After  acquiring  a  remote  terminal  and  a  tele¬ 
phone,  dial  the  system's  number  (FTS  217/333  4  086, 
WATS  800/637-0958),  If  there  is  no  answer,  the  entire 
system  is  down  for  maintenance.  Upon  hearing  a 
steady  tone,  plug  the  phone  into  the  terminal,  making 
sure  that  the  earpiece  and  the  speaker  are  in  the 
proper  openings,  and  log  into  the  system.  After  logging 
in  with  the  correct  name  and  password,  you  will 
receive  system  messages.  If  the  system  indicated  "NO 
DIRECTORY"  or  a  similar  designation,  access  to  ET1S 
is  probably  closed  down  and  you  should  try  again 
later.  If  the  system  is  in  operation,  output  similar  to 
that  shown  in  Figure  7  will  appear  on  the  screen  or 
terminal.  (The  text  on  the  right-hand  side  of  the  figure 
is  explanatory  and  is  not  part  of  the  output.) 

Keep  several  things  in  mind  when  using  the  system. 
The  symbol  <CR>  used  in  some  instructions  means  to 
depress  the  carriage  return  button.  The  instruction  to 
type  CTRL-d  means  to  simultaneously  depress  the 
button  marked  CTRL  and  the  letter  d.  An  input  error 
can  be  corrected  by  typing  CTRL-h  (hitting  the  CTRL 
and  h  buttons  simultaneously)  if  the  return  button  has 
not  yet  been  depressed.  This  procedure  will  back  the 
carriage  up  one  space  each  time  it  is  repeated.  This  can 
be  done  as  many  times  as  necessary.  Every  symbol 
which  has  been  backspaced  over  has  been  removed 
from  the  terminal  memory.  Therefore,  if  the  first 
digit  of  a  six-digit  number  has  been  mistyped,  you 
must  depress  CTRL-h  six  times  and  then  retype  all  six 
digits.  The  corrected  symbols  will  be  overprinted  on 
the  paper.  To  stop  a  long  listing  depress  the  button 
marked  DEL  (delete). 

The  Plant  Model 

Like  ET1S,  the  plant  model  also  contains  an  answer- 
driven  system  to  accept  required  input  data  from  the 
user.  Five  types  of  information  are  required  to  run  the 
model: 

1.  A  process  train  showing  all  flow  interconnec¬ 
tions,  source  water  inputs,  and  outfalls. 

2.  Volume  of  source  water  necessary  (cubic  meters/ 
hour  or  cubic  meters/batch  for  batch  processes)  at  a 
minimum  of  two  different  production  levels  (percent 
of  full  production). 


5S.  E.  Thomas,  R.  A.  Mitchell,  R.  E.  Riggins.  J.  J.  Fittipaldi, 
and  E.  W.  Novak,  Computer-Aided  Environmental  Impact 
Analysis  jor  Industrial,  Procurement,  and  Research,  Develop¬ 
ment,  Test,  and  Evaluation  Activities-  User  Manual,  Technical 
Report  N43/ADA056997  (CERL,  1978). 


3.  Volume  of  each  pollutant  (constituent)  gener¬ 
ated  (kg/hour  or  kg/batch)  in  the  water  phase  versus 
the  production  level  (percent  of  full  production). 

4.  Efficiency  (percent  removal)  of  all  treatment 
processes  for  each  pollutant  (constituent)  with  respect 
to  flow  rate  through  treatment  process  (cubic  meters/ 
hour  or  cubic  meters/batch). 

5.  How  the  process  train  could  be  modified  to 
incorporate  recycle/reuse. 

The  model  returns  values  for  the  flow  rate  and 
pollutant  loading  for  each  interconnection  in  the 
process  train,  including  the  outfalls.  Then  modifica¬ 
tion  of  the  existing  process  train  or  specification  of  a 
new  one  is  allowed  before  the  program  is  rerun. 

Several  general  instructions  should  be  followed 
throughout  the  program. 

1 .  All  process  names  are  alphanumeric  words 
which  can  contain  any  character,  but  should  not 
include  blanks  and  must  be  no  more  than  20  characters 
long.  It  is  important  to  maintain  spelling  and  upper/ 
lower-case  consistency.  It  is  also  important  to  make 
sure  that  a  process’  name  is  the  same  as  the  name  of 
the  process  model,  if  any,  that  is  on  file.  At  some,  but 
not  all  prompts,  the  program  will  find  a  spelling  error 
and  let  the  user  correct  it. 

2.  Yes  or  no  questions  can  be  answered  by  spelling 
out  yes,  no,  or  just  by  y  or  n.  Upper/lower  case  is  not 
important. 

3.  The  first  constituent  that  should  be  modeled  in 
any  run  is  water  because  water  flows  are  needed  by 
the  treatment  process  models.  Changes  in  flow  routing 
or  in  production  level  change  water  flow  rates,  and  the 
program  reminds  you  to  model  water  again  after 
making  these  changes. 

4.  At  any  point  in  the  program  where  a  list  is 
typed  in,  the  program  will  keep  repeating  the  same 
prompt.  Typing  a  0  (zero)  tells  the  program  that  the 
user  is  through  entering  the  list  and  wants  to  move  on 
to  the  next  prompt.  In  general,  “0"  acts  as  an  escape. 

5.  At  any  prompt  which  is  followed  by  “h/a," 
on-line  help  is  available.  Simply  answer  the  question 
with  the  word  “help"  to  obtain  a  help  routine  which 
provides  information  similar  to  that  given  in  this 
report. 
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ETIS  (Trademark  applied  for) 


United  States  Army  Corps  of  Engineers 
Environmental  Technical  Information  System 


ETIS:  What  program?  (Type  <cr>  to  see  list): 
Type: 


1  or  intro 


2  or  eics 

3  or  celds 

4  or  eil's 

5  or  afeics 

6  or  iicep 

7  or  xper 


8  or  help 

9  or  rubouts 
10  or  end  or  bye 

!  mail 


for  an  introduction  to  the  Environmental  Technical 

Information  System  and  people  to  call  for  help  with 
problems  about  using  the  programs, 
for  the  Environmental  Impact  Computer-aided  System, 
for  the  Computer-aided  Environmental  Legislative  Data  System, 
for  the  Economic  Impact  Forecasting  System, 
for  the  Air  Force  Environmental  Impact  Computer-aided  System, 
for  the  Interagency/Intergovernmental  Coordination  for 
Environmental  Planning  (IICEP). 
for  the  Experimental  Subsystem  Module  of  ETIS  - 
Pilot  Systems  under  development, 
for  HELP  on  using  any  of  the  ETIS  systems  or  UNIX, 
to  ignore  extraneous  phone  noise, 
to  exit  from  ETIS. 
to  see  your  mail. 


Figure  7.  Example  ETIS  instructions  and  explanation. 


Step-by-step  instructions  for  using  the  model  are 
given  below,  in  the  order  in  which  the  program  general¬ 
ly  flows.  Because  you  can  choose  how  the  program 
proceeds  at  several  points,  some  prompts  may  be 
repeated  during  a  session. 

1 .  Do  you  want  process  configuration  input  to 
come  from  a  file?  (h/a) 

This  allows  you  the  option  of  keeping  the  answers 
to  questions  3  through  5  (below)  in  a  file.  Be  sure  to 
include  all  the  zeroes  in  the  file  to  terminate  question¬ 
ing  (see  general  instruction  4  above). 

2.  What  is  the  name  of  the  file?  (h/a) 

Question  2  asks  the  name  of  the  input  file  if  ques¬ 
tion  1  was  answered  affirmatively. 

If  question  1  was  answered  negatively,  questions  3 
through  S  will  be  asked. 


3.  What  is  the  name  of  the  next  node?  (h/a) 

The  first  step  in  the  program  is  to  describe  the  flow 
network.  This  is  done  by  first  typing  the  names  of  each 
node,  in  any  order. 

4.  What  node  does  flow  go  to  from  <nodename>? 
(h/a) 

For  the  node  just  entered  in  response  to  question  3, 
enter  the  name(s)  of  the  node(s)  to  which  flow  goes. 
When  there  are  no  more,  type  in  "0." 

Questions  3  and  4  will  be  repeated  until  all  nodes 
are  entered;  at  this  point,  type  in  “0"  in  response  to 
prompt  3. 

5.  Enter  node  which  goes  to  sink,  (h/a) 

Type  in  the  name  of  any  process  from  which  flow 
goes  to  the  outfall.  This  prompt  will  be  repeated, 
allowing  multiple  outfalls,  so  type  “0”  after  the  last 
one. 
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Ai  this  point,  the  flow  connections  are  printed  out. 

6.  Do  you  want  to  change  flow  connections?  (h/a) 

To  change  where  flow  goes  for  any  of  the  nodes, 
answer  yes.  This  question  is  repeated  at  the  beginning 
of  each  run,  so  new  flow  routing  schemes  can  be  tried 
during  the  same  session. 

Prompts  7  through  9  occur  if  (low  connections  are 
changed. 

7.  Which  node  do  you  want  to  change?  Ih/a) 

Enter  the  name  of  the  node  whose  outflow  you 
wish  to  re-route.  This  instruction  will  be  repeated  at 
the  end  of  the  change  routine,  so  you  can  change  as 
many  nodes  as  desired. 

At  this  point,  the  nodes  which  flow  goes  to  from 
the  node  you  want  to  change  are  printed  out. 

8.  Which  node  do  you  wish  to  eliminate  flow  to? 
(h/a) 

If  you  want  to  stop  flow  from  going  to  a  node, 
enter  the  node's  name.  You  can  eliminate  flow  to  the 
outfall  by  entering  "outfall.”  This  question  will  be 
repeated  so  you  can  eliminate  flow  to  more  then  one 
node,  if  necessary. 

9.  Which  node  do  you  wish  to  add  flow  to?  (h/a ) 

If  you  want  to  direct  flow  to  a  node,  enter  the 
node’s  name.  If  not,  enter  “0.”  You  can  add  flow  to 
the  outfall  by  entering  "outfall."  This  question  is  also 
repeated. 

After  all  changes  have  been  made,  or  if  no  changes 
were  requested,  the  following  prompts  are  given. 

10.  What  production  level?  (h/a) 

Enter  the  production  level  you  want  to  model. 
This  should  be  expressed  as  a  percent  of  full  pro¬ 
duction,  so  it  should  be  between  0  and  100. 

1 1 .  What  constituent  do  you  want  to  model?  (h/a) 

At  this  point,  the  program  also  prints  out  a  list  of 
the  constituents  for  which  models  exist,  and  you 
should  choose  from  them.  You  may  also  add  a  new 
constituent  to  the  list  by  entering  its  name.  For  the 
first  run,  and  after  any  changes  in  flow  routing  or 


production  level,  water  should  be  modeled  first  The 
prompts  now  given  depend  on  whether  wa.er  or  a 
pollutant  is  being  modeled.  The  prompts  for  modeling 
water  flows  are  explained  in  questions  13  through  15. 
The  prompts  for  modeling  pollutant  generation  are 
explained  in  questions  19  through  28. 

1  2.  Do  sou  want  to  use  the  existing <consuluent> 
models'1  Ih  a) 

If  you  want  to  use  the  models  already  on  the  file, 
without  any  modifications  and  without  seeing  the 
models,  answer  "yes."  This  will  bypass  prompts 
13  through  16,  20  through  23.  and  25  through  28. 

13.  Is  there  a  water  demand  at  <m>denamc>'’ 
(h/a) 

If  the  process  represented  by  the  nodename  is  one 
which  requires  a  certain  amount  of  water,  answer 
"yes."  This  tells  the  program  that  you  wish  to  model 
the  water  demand  at  this  node. 

If  question  13  is  answered  Yes."  the  program  looks 
for  a  file  of  existing  models  which  includes  this  node. 
If  the  file  is  not  found,  or  it  no  model  with  the  node's 
name  is  found  in  the  file,  the  next  prompt  will  be 
question  15.  If  a  model  is  found,  it  will  be  graphed  out 
on  the  terminal. 

14.  Do  you  want  to  change  this  model? 

If  the  graph  adequately  represents  the  water  demand 
vs.  production  level  for  this  node,  answer  “no."  Other¬ 
wise.  answer  "yes,"  and  you  will  be  able  to  type  in  a 
new  model  by  answering  questions  1  5  and  16. 

15.  Production  level0 

16.  Water0 

Each  tin  e  these  two  prompts  are  given,  put  in  a 
point  on  the  water  demand  vs.  regular  production 
level  graph.  You  need  nui  enter  a  point  for  it ro  pro¬ 
duction.  When  a  production  level  of  100  is  entered, 
the  program  automatically  moves  on. 

The  model  is  now  graphed  ou  .  and  you  will  be 
asked  again  if  you  want  to  change  it.  When  you  are 
satisfied  with  the  model,  answer  "no"  and  the  new 
model  will  be  written  into  the  file. 

Questions  13  through  16  will  be  repeated  for  all 
nodes. 
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17.  Enter  fraction  of  flow  from  <nodename>  to 
<nodename>.  (h/a) 

Wherever  there  is  a  divided  outflow  from  a  node, 
the  program  must  know  what  fraction  of  the  flow  goes 
to  each  outflow,  ft  will  ask  this  question  for  all  nodes 
which  have  divided  outflows.  The  value  entered  should 
be  between  zero  and  one. 

After  question  17  has  been  answered,  the  program 
has  all  the  input  it  needs  and  will  compute  the  flows 
between  all  nodes.  The  results  will  be  typed  out  at 
the  terminal  as  well  as  appended  to  a  file  called 
RESULTS,  which  can  be  accessed  after  the  session  is 
over. 

1 8.  Do  you  want  another  run?  (h/a) 

Answering  “no”  will  terminate  the  session.  “Yes” 
will  take  you  back  to  question  6. 

The  prompts  explained  below  are  for  modeling 
pollutant  flows.  The  following  sequence  will  follow 
question  1 1  if  a  pollutant  constituent  is  chosen  for 
modeling. 

19.  Is  <constituent>  generated  at  <nodename>? 
(h/a) 

If  the  node  is  one  at  which  the  constituent  enters 
the  wastewater  stream,  answer  “yes,”  and  pollutant 
generation  will  be  modeled,  as  explained  in  questions 
20  through  23.  A  “no”  reply  will  lead  to  further 
prompts,  as  explained  in  questions  24  through  28. 
Replying  with  an  "x”  tells  the  program  that  this  node 
is  not  a  treatment  process;  this  bypasses  prompts  24 
through  28. 

After  a  “yes”  reply  to  question  19,  the  program 
searches  for  a  model  of  pollutant  generation  vs.  pro¬ 
duction  level  in  a  file.  If  no  file  exists  or  if  there  is  no 
model  in  the  file,  you  will  have  to  build  one  as  in 
questions  21  and  22.  If  a  file  does  exist,  it  will  be 
graphed  out  and  question  20  asked. 

20.  Do  you  want  to  change  this  model? 

If  you  are  satisfied  with  the  existing  model,  answer 
“no.”  Otherwise,  questions  21  and  22  will  be  used  to 
build  a  new  model. 

2 1 .  Production  level? 

22.  <constituent>? 


Put  in  the  points  for  a  graph  for  the  amount  of  the 
pollutant  constituent  generated  vs.  production  level. 
Production  level  should  be  expressed  as  a  percent  of 
full  production,  i.e.,  between  0  and  100.  You  need  not 
put  in  a  point  for  zero  production;  the  program  stops 
asking  for  points  after  you  have  entered  a  production 
level  of  1 00. 

23.  Do  you  want  to  change  this  model? 

If  you  are  satisfied  with  the  new  model  whicn  has 
just  been  graphed,  answer  “no,”  and  it  will  be  written 
into  the  file.  Otherwise,  you  will  go  back  to  question 
20. 

If  the  node  was  not  one  where  the  pollutant  was 
generated,  the  following  series  of  prompts  will  appear. 

24.  Is  <nodename>  a  treatment  process  where 
<constituent>  is  removed?  (h/a) 

Answer  “yes”  if  you  want  to  model  the  removal  of 
the  pollutant  at  this  node,  and  you  will  get  prompts 

25  through  28.  A  “no”  moves  you  on  to  the  next 
node. 

Modeling  treatment  processes  is  very  similar  to 
modeling  generation  or  water  demand,  except  you 
are  modeling  percent  removal  vs.  flow  rate  for  treat¬ 
ment  processes. 

At  this  point,  when  modeling  a  treatment  process, 
the  program  looks  for  a  file  containing  a  model  for  the 
process.  If  no  file  is  found  or  no  model  with  the  node's 
name  is  found  in  the  file,  questions  26  through  28  will 
be  used  to  build  a  model.  If  a  model  is  found,  it  will 
be  graphed  out  and  question  25  asked. 

25.  Do  you  want  to  change  this  model? 

If  you  are  satisfied  with  the  existing  model  of 
treatment  efficiency  (y-axis)  vs.  inflow  rate  (x-axis), 
answer  “no,”  and  the  program  will  use  the  model. 
Otherwise,  answer  “yes,”  and  you  will  get  questions 

26  through  28  to  build  a  new  model.  If  your  model 
does  not  include  flow  rates  as  high  as  you  are  trying 
to  model,  you  will  have  to  build  a  new  model. 

26.  Flow  rate? 

27.  Removal? 

For  up  to  20  points,  put  in  flow  rate  vs.  removal. 
Type  in  “0”  for  flow  rate  to  end  the  input  and  move 
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on.  The  new  model  will  be  graphed  for  you  at  the 
terminal. 


APPENDIX  A: 
EXAMPLE  SESSION 


28.  Do  you  want  to  change  this  model? 

If  you  are  satisfied  with  the  model  you  just  built, 
answer  “no,"  and  the  program  will  return  to  question 
19  until  all  nodes  have  been  examined.  A  "yes”  reply 
will  return  you  to  question  26. 

Appendix  A  provides  an  example  session  of  the 
model,  and  Appendix  B  provides  the  program  listing. 


4  SUMMARY  AND  RECOMMENDATIONS 

A  pilot  ammunition  plant  computer-aided  simu¬ 
lation  model  was  designed  and  programmed  to  predict 
wastewater  effluent  quality  at  various  production  levels 
and  to  assess  the  impact  of  reuse/recycle  on  effluent 
quality.  Any  type  of  process  train  can  be  used  as 
input  to  the  model  to  describe  water  use  in  AAPs.  The 
plant  model  may  be  accessed  through  the  Environ¬ 
mental  Technical  Information  System. 

Detailed  model  input  data  on  wastewater  quality 
and  quantity  for  cooling,  washing,  and  transport 
operations  is  not  available  in  the  literature.  Such 
information  is  probably  available  in  AAP  operating 
logs,  but  would  take  some  time  to  extract,  and  it  may 
not  be  available  for  more  than  one  production  level. 
Appropriate  estimations  and  assumptions  can  be  made 
by  the  user  when  water  use  and  pollutant  generation 
are  not  well-defined  functions  of  production  rate. 

It  is  recommended  that  the  plant  model  described 
here  be  tested  and  used  by  AAP  and  other  DARCOM 
personnel  before  being  further  developed  or  refined. 


This  example  run  is  for  a  simplified  version  of  the 
sellite  wash  process  for  TNT,  the  source  of  so-called 
“red  water."  For  this  example,  a  pond  is  the  only 
treatment  process.  In  recent  years,  “red  water"  has 
either  been  sold  for  paper  processing  or  dried  and 
incinerated.  In  the  example,  a  hypothetical  recycle 
of  the  pond  effluent  is  considered. 

The  model’s  level  of  resolution  can  be  chosen  by 
the  user;  nodes  may  be  chosen  for  units  as  small  as 
individual  reactors,  or  whole  process  trains  may  be 
included  in  one  “black  box”  node.  A  degree  of 
resolution  should  be  chosen  that  is  suitable  for  the 
data  available. 


sellite  mfgr- 


sellite  wash 


outfall 


Figure  Al.  Example  process  train. 
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Plant  Model  Program 

At  all  prompt*  which  are  followed  bg:  <h/a)  you  can  get  on-line 
help  bg  responding  with  the  word  'help' 

Answer  ges-no  questions  with  ges,  no.  g.  or  n 


Do  you  want  process  configuration  input  to  com*  from  a  file?  (h/ai 

Do  gou  want  general  help  or  help  specific  to  the  question  gou  are  on? 
Enter  1  or  2: 

1  =  general 
2=spec  if  1C 

© 


Ammunition  Plant  Process  Model 


This  is  a  program  to  modal  the  flow  of  water  and  water-borne  pollutants 
The  following  general  instructions  applg: 

1.  All  process  (node)  names  should  contain  20  or  fewer  characters  of 
ang  tgpe.  but  without  blanks.  Be  careful  With  spelling! 

2.  Ves  or  no  question*  can  be  answered  with:  ges.  no,  g,  or  n. 

3.  Whenever  a  prompt  repeats  itself  gou  can  tell  it  to  move  on 
bg  tgpinq  '0'  (zero). 

4  You  must  alwags  model  water  flows  before  pollutant  flow*  because  treatment 
models  use  flow  rat*  as  a  parameter.  Ang  time  gou  change  the  production 
level  or  re-rout*  flow*  gou  should  model  water  flows  again. 

5.  At  ang  prompt  where  '<h/a>'  appears,  gou  can  answer  'help'  ano  act  it. 

More  details  are  available  in  the  user  manual  and  details  are  available 
her*  for  specific  prompts. 


<8 


it  ' 1 '  to 


hoos*  one 


continue 

of  the  following: 


1  Do  gou  want  input  to  come  from  a  file? 

2  What  is  the  name  of  the  file? 

3  What  is  the  name  of  the  next  node? 

4  What  node  does  flow  go  to? 

5  Enter  node  which  goes  to  outfall 

6  Do  gou  want  to  change  flow  connections? 

7  Which  nod*  do  gou  want  to  change? 

8  Which  nod*  do  gou  wish  to  eliminate  flow  to? 

9  Which  nod*  do  gou  wish  to  add  flow  to? 

10  What  production  level? 

11  Uhat  constituent  do  gou  want  to  model? 

12  Do  gou  want  to  us*  the  existing  models? 

13  Is  there  a  water  oemand  at  this  node? 

15  for  help  with  units. 
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17  Enter  fraction  of  flow  from  a  to  b 
13  Do  you  want  anothar  run? 

19  I*  the  constituent  generated  hara? 

24  Is  this  a  treatment  process  where  constituent  is  removed^ 

98  for  intrapratation  of  results 

99  to  quit  HELP  and  return  to  program 

Which  one  do  you  want?  (you  entered  'help'  from  no  1 )  (T) 

Help  for  question  1: 

Do  you  want  process  conf igurati on  input  to  come  from  a  file? 


If  you  hove  a  file  with  all  tha  node  names,  where  flow 
goes  to  for  aach  one.  and  the  names  of  nodes  which  dram 
to  the  outfall,  you  can  use  it  instead  of  answering 
questions  3  through  5.  This  is  vary  useful  whan  running  a 
model  several  times.  If  you  answer  yes,  tha  neit  question 
will  ask  you  tha  name  of  the  input  file 

More  help?  y  or  n  (r^ 

Do  you  want  process  conf iguration  input  to  coma  from  a  file?  <h/a>(n) 

Process  configuration  entry  routine,  put  in  the  name 
and  where  flow  goes  to  for  each  node 


sell itemfr 
What  node  does 
sell i tewash 
What  node  does 

0 

What  is  the  ns 
sc  1 1 1 tewash 
What  node  does 
treatmentpond 
What  node  does 
0 

What  is  the  ns 
treatmentpond 


What 

is  th 

e  name 

of  th 

•  n 

est  node? 

0 

Ent*r 

node 

which 

goes 

to 

sink  <h/a> 

Enttr 

node 

which 

g09  9 

to 

sink  (h/a) 

Plow  goes  from  sellitemfr  to  sellitewesh 
Plow  goes  from  sellitewesh  to  treatmentpond 


Flow  goes  from  treatmentpond  to  sink 

Do  you  want  to  change  flow  connections?  y  or  n  (h/a)  (n) 
What  production  level-  0  to  100?  <h/a)  ^Qo) 

Whet  constituent  do  you  want  to  model? 

The  following  are  available:  <h/a) 
water 
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Do  you  want  to  use  the  existing  water  models'*  <h/a)@ 
Is  there  a  water  demand  at  sal  1 itemf r?  y  or  n  <h/e>  Q) 
Model  Water  Use  at  sellitemfr 
No  file  of  models  found;  file  created 


BUILD  A  NEW  MODEL  for  sell itemf r 

Input  water  used  as  a  function  of  production  level 

ror  a  max  of  20  points,  put  in  prod  level,  up  to  ICO, 
not  including  zero,  and  water  used/ senerated 

Use  units  of  cubic  meters/hour  or  cubic  meters/batch  for  water 
(for  help  with  units  type  'help') 

P^od.  l*velCl3?  (h/a)  (f5) 

water?  (2^ 

Prod  l#v*U2]?  (h/a)  @) 
water?  (^5) 

Prod  leve  1  C3]'>  <h/a>  (lOo) 


water"*  (105 
105  00 


64  00 


63  00 


42  00 


21  00 


60 


90 


20  40 

PRODUCTION  LEVEL 
Do  you  want  to  change  this  model"?  y  or  n  Q 

Is  there  a  water  demand  at  sel  1  i tewesh?  y  or  n  (h/a>  (n) 

Is  there  a  water  demand  at  treatmentpond?  y  or  n  (h/a)  (n) 


100 


water  from  sellitemfr  to  sellitewash  ■  105  00 


water  fr am  sellitewash  to  treatmentpond  •  105  00 

water  from  treatmentpond  to  outfall  »  105  00 

Inflow  to  network  at  sellitemfr  *  105  00 


Do  you  want  another  run?  y  or  n  (h/a)^ 


Do  you  want  to  change  flow  connections?  y  or  n  (h/a)  (n) 

What  production  level-  0  to  ISO'!'  (h/a) 

(Previous  level  was  100)  (iOu\ 


what  constituent  do  you  want  to  model? 
The  following  are  available,  (h/a) 


DNT  (TNT 


Do  you  want  to  use  the  existing  TNT  models?  (h/a)^^ 


Is  TNT  generated  at  sellitemfr?  y  or  n  (h/a)  (he^t| 

Do  you  want  general  help  or  help  specific  to  the  question  you 
Fnter  1  or  2: 


1«  general 
2«specif ic 

© 

Help  for  question  19.  Is  constituent  generated  here 
If  this  is  a  process  at  which  the  pollutant  enters  the 
wastewater  stream,  and  at  which  there  is  a  pollutant 
generation  vs.  production  level  model,  answer  yes 
Answering  with  a  tells  the  program  that  there  is 

neither  a  pollucant  generation  nor  a  treatment  model 
for  this  node. 

More  help?  y  or  n  (n) 

Is  TNT  generated  at  sellitemfr?  y  or  n  (h/a)  Q 
Is  sellitemfr  a  treatment  process  where  TNT  is  removed?  (h/a) 
Is  TNT  generated  at  sellitewash?  y  or  n  (h/a)  ^ 

No  file  of  models  found;  file  created 


BUILD  A  NEW  MODEL,  for  sellitewash 

Input  TNT  generated  as  a  function  of  production  level 

For  a  max  of  20  points,  put  in  prod,  level,  up  to  ICO, 
not  including  zero,  and  TNT  used/generated 

Use  units  of  kg/hour  or  kg/batch  far  TNT 
(for  help  with  units  type  'help') 
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Prod  ltv«lCl3?  <h/*> 

TNT''  (0) 

Prod  lava 1 C23?  <h/a)  (0) 
TNT"’  200 

Prod  levelC3P  (ft/a)  1100 


TNT'  (250) 

Prod  l»v*lC4]?  <h/a>  (100) 
TNT?  (loo) 

250.  00 


200  00 


150  00 


100.  00 


50  00 


oo 

000  o 
00  o 

000  0 
0  0 
00  0 


oooooo 


-I  — 

to 


20  40 

PRODUCTION  LEVEL 

Do  you  want  to  changa  Uit  modal?  y  or  n  ^ 


SO 


100 


BUILD  A  NEW  MODEL  for  sallitaieash 

Input  TNT  ganaratad  at  a  function  of  production  laval 

For  a  max  of  20  points*  put  in  prod.  laval*  up  to  100* 
not  mcludtng  taro,  and  TNT  usad/ganaratad 

Use  units  of  kg/hour  or  kg/batch  for  TNT 
(for  help  (pith  units  typa  'halp') 


Prod.  tavalCU?  (h/»>  (10) 
TNT?  0) 

Prod.  lavaIC2J?  <h/a>  0) 

TNT">  (200} 

Prod,  laval C3J?  (h/a>  0 
TNT">  (0 

Prod.  Iavalt4j?  <h/a)<J3£> 
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TNT?  (30^ 
300  00 

240  00 

130. 00 

120  00 


60  00 


ooo 


60 


80 


100 


20  40 

PRODUCTION  LEVEL 

Do  you  want  to  change  this  model?  14  or  n  (r^ 

Is  TNT  genera ted  at  treatmentpond?  y  or  n  th/a)  @ 

Is  treatmentpond  a  treatment  process  where  TNT  is  removed?  lft/a> 

Mo tel  Treatment  Efficiency  at  treatmentpond 
No  model  called  treatmentpond  in  file 

BUKO  A  NEW  MODEL  for  treatmentpond 

Put  in  removal  percent  as  a  function  of  flow 

For  a  ma»  of  20  points,  put  in  flow  rate  and  removal 

Use  units  of  cubic  roeters/hour  or  cubic  meters/batch  for  flew 
For  help  with  units  type  'help' 

Flow  rateC  1 17  <h/e>  (30) 

removal?  ^8) 

Flow  rate£2I?  (h/#>  (kw) 

removal? 

Flow  rateC33?  <h/a>  (5O0) 
removal?  (75) 

Flow  rateC43?  th/a)  (5) 


T 

R 

E 

A 

T 

M 

E 

N 

T 


<98.  00 


78  40 


58  80 


00000900 


0090000S9C 


0000000000 


0000000000 
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39  20 


i 

: 


r 

e 

N 

C 

V 


19  60 


100  200 
FLOW  RATE 

Do  you  want  to  change  this  model"7  y  or 


300 

"  © 


400 


500 


TNT  from  sellitemfr  to  sellitewash  a  0.00 
Concentration  =  0  mg/1 

TNT  from  sellitewash  to  tr eatmentp ond  »  300.00 

Concentration  a  2857  14  mg/1 

TNT  from  treatmentpond  to  outfall  =  16  00 

Concentration  *  171  429  mg/1 


Do  you  uiant  another  run?  y  or  n 


<  h/a  >  © 


Do  you  want  to  change  flow  connections?  y  or  n  (h/a)  (7 


Wh  ch  node  do  you  want  to  change?  (h/a) 

Do  you  want  general  help  or  help  soecific  to  the  question  gnu  are  or, 
Enter  1  or  2: 

1=  general 
2=spec l f  i  c 

© 

Help  for  question  7:  Which  node  do  you  want  to  change? 

You  are  now  in  the  routine  to  change  flow  routing 
Enter  the  name  of  the  node  (process)  whose  outflow  you 
wish  to  red  i  rec  t 

More  help7  y  or  n 


Which  node  do  you  want  to  change?  (h/a)  treatment 

No  such  node  as  pond 

Which  node  do  you  want  to  change?  (h/a)  (jreatmentponi}) 
Flow  goes  from  treatmentpond  to  outfall 
Which  node  do  you  wish  to  eliminate  flow  to?  (h/a)© 


CjonT) 


Which  node  do  you  wish  to  add  flow  to?  (h/a)( 


don? 
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Which  nod*  do  you  with  to  add  flow  to'7  <h/a>(^ 


Which  node  do  you  want  to  change7  (h/a)(^ 

Plows  changed-  model  water  first  to  get  correct  results 
What  production  level-  0  to 100?  (h/a) 

(Previous  level  was  100)  Cl 00} 

What  constituent  do  you  want  to  model? 

The  following  are  available:  (h/a) 

water 


TNT 

ON  T  ^uater) 

Do  you  want  to  use  the  existing  water  models?  (h/a)(u) 

Is  there  a  water  demand  at  sellitemfr"*  y  or  n  (h/a)  Q 

Is  there  a  water  demand  at  sellltewash?  y  or  n  (h/a)  (n) 

Is  there  a  water  demand  at  treatmentpond?  y  or  n  (h/a)  (ri) 

Enter  fraction  of  flow  from  treatmentpond  to  outfall ^8) 


water  from  treatmentpond  to  sellitemfr  »  22.00 

water  from  sellitemfr  to  sellltewash  «  i  !•  00 

water  from  sellitewash  to  treatmentpond  *  1 lO.  00 

water  from  treatmentpond  to  outfall  *  88.00 

Inflow  to  network  at  sellitemfr  »  98.00 


Do  you  want  another  run?  y  or  n  (h/a)^ 


Do  you  want  to  change  flow  connections?  y  or  n  (h/a)  Q 

Whet  production  level-  0  to  ipO?  (h/a) 

(Previous  level  was  100)  GocQ 

What  constituent  do  you  want  to  model? 

The  following  arc  available  (h/a) 

watar 


TNT 


Do  you  want  to  usa  tha  existing  TNT  models?  <h/a>(^ 
Is  TNT  generated  at  sellitemfr?  y  or  n  (h/a)  Q 
Is  TNT  generated  at  sailitewash?  y  or  n  (h/a)  ^ 
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I 


Is  ~NT  generated  at  trea tm*n tp ond  “  y  or  n  (h/a>  © 

Is  treatmentpond  a  treatment  process  where  TNT  is  removed?  (h/'a)Q 
Enter  fraction  of  flow  front  treatmentpond  to  outfall 


TNT  from  treatmentpond  to  sellitemfr  *  3  64 

Concentration  =  165  624  mg/1 

TNT  from  sellitemfr  to  sellitewash  =  3.64 

Concentration  =  33. 1243  mg/1 

TNT  from  sellitewash  to  treatmentpond  *  303  64 

Concentration  —  2760  4  mg/1 

TNT  from  treatmentpond  to  outfall  =  14  57 

Concentration  =  165  624  mg/1 


Co  you  want  another  run?  g  or  n 


<  h  /a  )  ^ 


Do  gou  want  to  change  fKo  connections'1  g  or  n  (h/a)  Q 

4hat  production  level-  0  to^iOO'  (h/a) 

(Previous  level  was  100>  (50) 

Production  level  changed-  model  water  first  to  get  correct  results 
Uhat  constituent  do  gou  want  to  model? 

The  following  are  available  (h/a) 

water 


TNT 

DNT  (tuaTer) 

Do  gou  want  to  use  the  enisting  water  models?  (h/a>(u) 

Is  there  a  water  demand  at  sellitemfr?  g  or  n  (h/a)  @ 

Is  there  a  water  demand  at  sellitewash?  g  or  n  (h/a)  (rt) 

Is  there  a  water  demand  at  tr eatmentpond"*  g  or  n  (h/a)  @ 

Enter  fraction  of  flow  from  treatmentpond  to  outfall 


water  from  treatmentpond  to  seilitemfr  »  13.00 

water  from  sellitemfr  to  sellitewash  *  65.00 

water  from  sellitewash  to  treatmentpond  ■  65.00 

water  from  treatmentpond  to  outfall  »  52.00 

Inflow  to  network  at  sellitemfr  =  52.00 


Do  gou  want  another  run?  g  or  n  (h/a>^ 

Do  gou  want  to  change  flow  connections?  g  or  n  (h/a) 
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What  oroduction  lev»l-  0  t 
i  Previous  lavtl  was  50)  ( 


100">  (h/a) 


What  constituent  do  you  want  to  model? 
The  following  are  available:  (h/a) 
water 
TNT 


Do  you  want  to  use  the  existing  TNT  models^  <h/a)Q 
Is  TNT  generated  at  sellitemfr?  y  or  n  (h/a) 

Is  TNT  generated  at  seilitewash?  y  or  n  (h/a)  © 

Is  TNT  generated  at  tr  ea  tmen  tp  ond"7  y  or  n  (h/a)  Q 


Is  treatmentpond  a  treatment  process  where  TNT  is  removed"7  (h/a) 


Enter  fraction  of  flow  from  treatmentpond  to 


outfall  Q 3 ) 


TNT  from  treatmentpond  to  sellitemfr  =  1  01 

Concentration  »  78  0065  mg/1 

TNT  from  sellitemfr  to  sellitewash  =  1.01 

Concentration  =  15  6013  mg/1 

TNT  from  sellitewash  to  treatmentpond  »  169  01 

Cone antrat i on  =  2600  22  mg/1 

TNT  from  treatmentpond  to  outfall  =  4.06 

Concentration  *  70  0065  mg/1 


Uo  you  want  another  run?  y  or  n  (h/a)^^ 
12  »  W 
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5  jj  C  (caTfiESULT^) 

next  run 

Constituent  water  Production  level.  100 
water  Prom  sellitemPr  to  sellitewash  =  105  CO 

water  Prom  sellitewash  to  treatmentpond  =  105  00 

water  from  treatmentpond  to  outPall  =  105  00 

InPlow  to  network  at  sellitemPr  =  105.00. 

NEXT  RUN 

Constituent  TNT  Production  level:  100 

TNT  Prom  sellitemPr  to  sellitewash  =  0  00 

Concentration  =  0  mg/1 

TNT  Prom  sellitewash  to  treatmentpond  >»  300  00 

Concentration  =  2857  14  mg/1 

TNI  Prom  treatmentpond  to  outfall  =  18  00 

Concentration  »  171.429  mg/1. 

NEXT  RUN 

Constituent  water  Production  level  100 
water  Prom  treatmentpond  to  sellitamfr  •>  22  CO 

water  Prom  sellitemPr  to  sellitewash  *  110.00 

water  Prom  sellitewash  to  treatmentpond  =  110.00 

water  Prom  treatmentpond  to  outPall  =  88  00 

InPiow  to  network  at  sellitemPr  *  88.00 

NE*T  RUN 

Constituent  TNT  Production  level  100 

TNT  Prom  treatmentpond  to  sellitemPr  =  3  64 

Concentration  «  165  624  mg/1 

TNT  Prom  sellitemPr  to  sellitewash  «  3.64 

Con  :  er.trut  ion  »  33  1248  mg/1 


TNT  Prom  sellitewash  to  treatmentpond  »  303  64 

Concentration  »  2760  4  mg/1 

TNT  Prom  treatmentpond  to  outPall  *  14  57 

Cone entrat i on  »  165  624  mg/1. 

NEXT  RUN 

Constituent  water  Production  level  50 
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13  00 


u»at*r  from  treatmentpond  to  sellitemfr  * 
water  from  sellitemfr  to  sellitewash  «  65  00 

water  from  sellitewash  to  treatment? ond  =>  65  00 

water  from  treatmentpond  to  outfall  »  52  00 

Inflow  to  network  at  sellitemfr  »  52  00 

NEXT  RUN 

Constituent  TNT  Production  level  50 

TNT  from  treatmentpond  to  sellitemfr  »  1  01 

Concentration  =  78  0065  mg/1 

TNT  from  sellitemfr  to  sellitewash  =  1  01 

Concentration  »  15.6013  mg/1 

TNT  from  sellitewash  to  treatmentpond  »  169  01 

Concentration  »  2600  22  mg/1 

TNT  from  treatmentpond  to  outfall  =■  4.06 

Concentration  »  78.0065  mg/1. 
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APPENDIX  B: 
PROGRAM  LISTING 


•♦include  ;stdio  n>; 

double  determCl 00 3  C 1 00 3 , 
double  r  h  s  C 1003 , 
douole  var C 1003 . 
int  d  e  trms  2 . 

int  matr  i  x  C303  C303.  argsiti 
struct  a  {char  aaC203,  }, 
struct  {char  ndnameC2G3. 

struct  a  ndtoC63;  ' 

double  wa  t  r i n i 
nodeC303i 

ma  in  <  ; 

{ 

int  nvar,  i.  j.  varlblC1003i 

int  li.  lx.  temp.  col.  row.  flag. 

int  noden,  ncnsts.  fdi 

int  tonodeC63i 

double  f 1 o  f rac  C63 .  conci 

double  qreq.  rmval,  qtrti 

double  flowsC1003, 

float  pi.  ftemp; 

char  tqueryC43.  queryC53.  ucnsttC203.  linebufC813; 
char  stempC53i 

struct  {char  aC203.  >  cnstsC203; 

/*  */ 

'*  Bildsys  is  the  main  routine.  It  calls  routine  MatTxin 
to  build  a  connection  matrix  for  the  system 
Routine  chgmtrx  allows  the  user  to  revise  the  c  or.r.ec  ti  or. 
matrix  Routine  findmdl  uses  a  file  of  production  level 
vs  water  demand  models  to  determine  how  much  water 
is  needed  at  each  nodeiplant  process)  which  requires  water 
Bildsys  uses  the  above  information  to  build  a  system 
of  simultaneous  linear  equations  describing  flow  through 
the  system.  Additional  equations  are  added  for  each 
node  in  the  system  with  more  than  one  outflow. 

After  the  system  is  built,  with  the  equations  described  bu  the. 
determinant  "determ"  and  right-side  vector  ”rhs".  it  is 
salved  in  the  routine  "solvsys 


*/ 

/'*  variables 
/  *  nvar 

/* 

/* 

/*  arysze 

/* 

/*  matrix 

determ 
rh  s 

var 1 b  1 

/*  ***a«4««««4ax#****»ax»*****a*t-'A**4*x«*t*t4  4«i  */ 


/*  print  out  heading  »/ 

pr 1 n t f ( " \n\n\n  Plant  Model  Program  ); 

printf<"\n\n  At  all  prompt*  which  are  followed  by  (h/a)  you  can  gee  on-line 
p r i n t f < " \nhe 1 p  by  responding  with  the  word  'help'  \n\nAn<wer  yet  no  ">• 
pr  i  ntf  (  "quest  i  on*  with  ye*,  no.  y.  or  n  \n\n\n\n">; 

/* _ _ _ __ _ _____ _ _ _ _  */ 


*/ 

the  number  of  simultaneous  equations  */ 
originally  equal  to  the  number  of  */ 
positive  values  in  the  connection  matrix  */ 
the  length  of  the  x(“from")  dimension  in  the  * 
connection  matrix  */ 

the  connection  matrix 

the  determinant  of  the  system  of  simultaneous  equation* 
the  right  side  (constants)  for  the  equations 
labels  variables  in  the  determinant  •/ 


/*  initialize  list  of  constituents 


r  -- 


*/ 


setbuf(stdout. NULL ) ; 
strc  py  (  cnstsCOl.  a.  “water"); 
strcpytcnstsCll  a. "TNT"); 
strcpy(cnstsC23  a. " DNT ” > ; 
ncnsts  =  3. 


/*  <<<<<<{<<<<<<<<<  Initialize  watrin  >>>>>>>>>>>>'.>>>>>>>>>  <<■ 
for  (i  =0;  i  <»29;  i++;  node!  i  3  watrin  *  0.0; 

P 1  *  -1 . 

/'*  /)));'/>;»»  call  input  routine  ««<^< '««<  "  :««  */ 

arysze  *  matr  i  in  <  ) , 

■'*  :  :  >:  :•:•»»>»»  print  and  modify  matrix  <<<<<<<<C<<<  <«««  */ 

for  (col  *  1,  col  <*  arysze;  col++) 

< 

for  (row  ■  1;  row  <■  arysze;  row**) 

{  i  f  (matr  i  xtrowK col  3  !*  0) 

pr  intf  (  "\nFlow  goe*  from  %s  to  */.s".  nodeCcoi  j  ndname. 
nodeCrowl.  ndname); 
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pr  i  nt  f  ( 11  \n 11 )  , 


for  (col  =1.  col  <=  arysze.  col++) 

if  wnatr  i 1  Car  y  s  ze+1  He  ol  3  !=  0) 

p r  i  n t f < " \nF 1 ow  goes  from  Xs  to  si  nk  " .  nod  e  C  c  o  1 1.  ndname  >. 


recycle  pr intf < " \n\nDo  you  want  to  change  flow  connections"  u  or  n  (h/a)  "  )  - 

sc  an  f  < "Xs " . quer y )  ; 

if  ( stremp ( query,  "help " )  ==0  !!  stremp (query.  "HELP" )  ==  0  ) 

(  help(t)i 

goto  recycle; 

> 

if  (querycOl  ==  'y'  i!  aueryCO]  *«*  'Y') 

<  c  hgmt  > ( ) ; 

pr intf ( "\nFl out  cnanged-  model  water  first  to  aet  correct  results”); 

> 


ten  pr intf < "\nWhat  production  level-  0  to  100?  (h/a)  ">; 
if  (pi  >  0  0)  printf  (  "\n(Previous  level  was  X-.  Of)  ",  pi); 
scanf  (  "V.s  "  ,  query  ) ; 

if  ( stremp ( query. "help" )  **  0  ! ;  stremp ( query, "HELP" )  “0 
<  help(lO); 
goto  ten; 

> 

sscanfiquery, "Xf " . if t emp ) , 

while  (ftemp  <0.0  !!  ftemp  >  100  0) 

<  pr intf ( "NnPlease  keep  production  level  between  0  and  100"), 
printf < "\n\nProduction  level?  ">; 
scenf  (  ”Xf  &f temp  ) ; 

> 


/*  if  the  production  level  is  being  changed,  say  so  */ 
if  (ftemp  pi  &&  pi  -10) 

pr i ntf (" \nProduct ion  level  changed-  model  water  first  to  get  correct  results"); 
pi  =  ftemp; 


i  =  99. 

while  (i  >«  nensts)  /*  read  in  constituent  wanted  and 

make  sure  it  is  available  •/ 

{ 

eleven  pr  int  f  (  "  VnWhat  constituent  do  you  want  to  model?"), 
printf ( ”\n\nThe  following  are  available:  (h/a)  ”>; 
for  < i»0>  i  <  nensts;  i++> 

<  printf  <  "\n\nXs  " ,  cnstsC  i  3.  a  ) ; 
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> 

sc  an  f  (  "Y.*  "  ■ ucnstt); 

if  i  s  trcmp  <  ucnstt.  "help  “  >  =-  0  strcmp (ucnstt. "HELP" >  == 

{  heip  <  1 1  ) . 
goto  eleven; 

> 

for  (i=0,  l  <  ncnsts;  i++) 

<  if  (  strcmp  <  cnstsC  i  3.  a,  ucnstt )  ==  0)  break; 

y 

if  (  1>=  ncnsts)  pr  intf  (  "\n'/.s  not  avai  1  ab  1  e  " ,  ucn  s  1 1 1 , 


,  v.  v.  y.v.  v.  y.y.  y.  y.  v.  v. 'L7.  y.y.  ■/.■/. ‘i'/.'/./.a'/. 


Allow  the  user  to  uce  the  existing  process  models  without 
mod i f i cat i on 

twelve  printf("\nDo  gou  want  to  use  the  existing  '/.%  models'?"  (h/a)  ",  ucnstt 
scanf ( "Xs” , tquer  g  ) . 

if  <tquery[03  ==  'Y')  tqueryiG3  «  'y'i 

if  ( strcmp ( tquerg,  "HELP" )  ==  0  !  !  strcmp ( tquerg,  "help  "  >  ==  0  ) 

<  help  ( 12); 
goto  twelve; 

> 


/*  if  modeling  water,  reset  the  stored  values  of  water  flows  to  zero 

i  f  <  s  trcmp  (  ucnstt,  "water "  )  ==  0)  <  for  (i=0;  i  O  29;  i++) 

-.:deC  i  3.  watr  in  =  0  0, 


> 


*/ 


All  the  input  and  bookkeeping  is  done  Now  smld 
the  system  of  simultaneous  equation 


////////// 


nvar=0> 

for  <i  =  l,  iO<  arysze+1 ) ;  i +••*•>  < 
for(j=l,  jOarysze;  j*+)  < 

l  f  <  ma  tr  i  x  C  i  3  C  j  1  !=0>  matrixti3Cj3  »  renvan 
>; 


*  array  now  read  into  "matrix"  */ 

/*  now  generate  determinant  for  system  of  equations  */ 
/*  put  "1"  in  determ  for  inflows  */ 

/*  and  "-1"  in  for  outflows  */ 

for  <i=0,  i  ‘>99;  1++) 
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\  for  i  j=0;  j  ^'=  99,  j++)  determi 1 3C j 3  =  0  0, 
rhsC  l  3  =  0  0; 


for  i  i=l,  i<=  arysze;  i++)  < 

for  ij  =  l,  j<=  arysze,  j++)  < 

if<  iMtrn  CilCjl  '=0  )  { 

d  e  ter  mC  i  3  C  ( ma  tr  i  j  C  i  3  f  i3)3  -  1  0, 
d etermC j 3C < matr l « C i j [  i  3 > 3  -  -1  0. 

>,  >,  >, 

for  <  j  =  1 ,  j<=  arysze,  j++  ) 

{ 

if(  matri«C  <a»ysze+l)  3Cj3  !=  0  ) 

determC  j  DC  (matr  i  x  C  <arysze+l>  3tj3>  3  =  -1.0; 

> 

decrmsz  =  arysze, 

/»  add  equations  for  inflows  to  siistem-  */ 

.  «•  one  equation  for  each  node  with  a  production-level  vs  deir.ano  function  #/ 
*  for  water  */ 

»  call  routine  "findmdl"  to  get  information  on  these  nodes  */ 

/'*  qreq  is  the  flow  rate  demanded  at  the  node  •/ 

.*  add  a  new  row  in  tne  determinant  for  each  inflow  */ 

fori  1  =  1,  i  .>argsie,  1-rr) 

< 

if  (  strcmp  (ucnstt.  "water"  )  **•  w> 

{ 

ql3  query C03  ■  'n'; 

printfi  "\nls  there  a  water  demand  at  7.s?  y  or  n  (h/a)  ",nod*Ci3  ndname). 
scanfi  "%s  ",  query); 

if  ( strcmp ( query, "help" )  •=  0  !!  s trc mp ( query , "HELP" )  == 

<  help!  13), 

goto  ql3; 

> 


i f  v  query  C03  ==  'Y'  !1  query[03  ==  'y') 


i *  build  new  equation  if  there  is  inflow  */ 

< 

ifitqueryt03  !=  'y  '  ) 

pr  int  f  <  "\nMod  ol  Water  Use  at  V.sXn" ,  node!  i  3  ndname  ) , 

/•  call  findmdl  with  a  1  to  allow  changes, 
with  a  zero  for  no  changes  •/ 

ifitqueryCQl  '»  'y  '  > 

qreq  ■  findmdl  <  nodeC  i  3.  ndname,  ucnstt,  p  1,  1 ) ; 
else  qreq  =  f indmd 1 ( nod eC i 3 .  ndname, ucnstt ,  p 1 . 0 ) ; 

++detrmsz; 

*+nv»r; 

varlblinvarl  =  i; 

for  <  j=l;  j<=arysze,  j++  ) 
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{  if(  uatrutiKjl  !=  0  ) 

determl det rms 2 3 C < matr 1 x C 1 3 t j 3 )  3=1.  >. 

d a t ermCde trmsz 3  Invar  3  =  1.. 
rnsCdetrmsz3  =  qreq, 

find  the  first  vsr  in  the  "ude's  column  */ 

ii  =  l. 

while  (  matrixCil3Ci3  =—  0  )  11++; 
for  (  ii=l,  ix<=detrmsz;  ix++  ) 

{  if<  determCixlC  <  ma  tr  i  x  C  i  i  3  C  i  3  )  3  ==  —1  0  >  determt  i  x  3  Cnvar  3  =  1  u;  >. 

> 

>  /*  end  of  water  inflow  routine 

Now  modify  equations  for  nodes  at  which  pollutants 
are  generated  Change  equations  from 
"Sum  of  inflows  ==  Sum  of  outflows”  to 

"Sum  of  inflows  ♦  pollutant  generated  ==  sum  of  outflow 
This  is  done  by  putting  -1  *  pollut  generated  on  right 
side  of  equation  *' 


else 

{ 

ql9  query  C03  =  'o'; 

pr  int  f  (  “ \nls  Xs  generated  at  Xs?  y  or  n  (h/a)  ,  ucnstt.  node!  i  3  ndname) 

scanf  (  "Xs",  query); 

i f ( strcmp ( query, "help " )==0  !1  str cmp ( query . "HELP” >  *-  0) 

{  he  1 p ( 19 ) ; 
goto  ql9, 

> 

if<queryC03  ==  'V'  !!  queryC03  ==  'y'  ) 

{ 

l f  ( tquery [03  ! =  '  y ' ) 

qreq  =  f indmd 1 < nod eC i 3  ndname. ucnstt , p 1 , 1 ) ; 
else  qreq  =  f i ndmd 1 < nod eC i 3  ndname. ucnstt. o 1 . 0 ) ; 
rhsCi)  =  -1.0  *  qreq; 

> 


•»  Modify  those  equations  for  noaes  which  are  waste  treatment  processes 
By  multiplying  the  "-l”  values  in  the  equation  by  the  fractional 
removal  by  the  treatment  process,  the  equation  is  changed  f’-em 
"inflow  ==  outflow"  to  "inflow  *  removal  *=  outflow" 

This  is  done  only  for  pollutants,  not  for  water  */ 


else 

< 


if  ( query  C  03 


if  the  user  replied  to  the  previous  ouestion 
with  a  assume  this  is  not  a 

treatment  process 

'x'  $<8<  queruC03  '=  'X'> 


37 


<  q24  p  r  int  f  <  "  \n  I  s  5Cs  a  treatment  process  where  /.s  is  removed""  <h/'a> 
nodeC i 3  ndname, ocnstt), 

sc anf  (  '’Xs",  queru ) , 

if  ( strcmp ( query, “help” )==0  s tremp ( query ,  "HELP " >  =~  0 
{  h  e 1 p 1 24 ) ; 
goto  q24; 

> 

> 

if  iqueryCOl  ==  'Y'  !!  queryCOi  =='y'  ) 

< 

if  (tqueryC03  1 =  'y') 

pr  intf  (  "\nModel  Treatment  Erficiency  at  ‘/.s  " ,  nod  eC  l  3  ndna/ne). 

/*  make  sure  water  has  been  modeled  first  and  theTe 

is  a  water  inflow  to  this  treatment  nod#  know/"  */ 

if  (nodeC  i  3 .  wati  in  ==  0.01 

{  pr intf < "VnWarning-  Water  flow  to  process  %s  is  Uro-  ", 
nodeC i 3  ndname ) i 

printf < "NnWere  water  flows  modeled  first??”), 

> 

y*  if  changes  in  model  are  allowed,  call  trtmdl  with  a  1 
if  not  ■  call  it  with  a  zero  */ 

i f  <  tquery  CO]  1 =  ' y ' ) 

rmval  =  tr  tmd  1  < nodeC  i  3.  ndname.  gens  1 1,  nodeC  l  3  watr  i n,  l  ) , 
else  rmval  =  trtmd  1  ( nodeC  i  3  ndname,  uenstt,  nodeC  l  3.  watr  in,  0,  • 

rmval  =10-  (rmval  /  100  0); 
for  (col=l,  col  O  nvar,  col++) 

C 

if  <  d etermC  i  3Ccol  3  >0.0  )  de t ermC l  3 C c o  1  3  » 

determC i 3Ccol 3  »  rmval. 

} 

> 

> 

> 

> 

* 

add  equations  for  each  node  with  divided  outflows 
telling  how  much  flow  goes  out  each  path, 
noden  is  the  node  number  which  has  divided  outflows 
"tonode"  is  a  matrix  of  the  nodes  flow  goes  to 

"flofrac  is  a  matrix  of  the  fraction  of  the  outflow  to  the  node 

/ 


for<col  =  1.  col  v=  arysie,  col+e) 

<  flag  =  0, 

for  (row  =1,  row  <*  arysxerl,  rower) 
{  l f ( matr i x Crow3 C col  3  1 =  0) 

{  if  ( Flag  1 =  0) 

<  noden  =  col. 

tonodetflaj3  =  ? ow. 
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ql7  if(ruw  ==  argsie  +  1  ) 

pr  in  t  f  <  “ \r.Ent  er  fraction  of  flow  from  to  outfall 
nodetcoll  ndname). 
rise 

pr  i  n  t  f  (  “\nEn  t  er  fraction  of  flou  from  /.s  to  v.s  <h/a> 
nodeCcoii.  ndname.  nodetrowl.  ndname); 
scanf  ( ”7.s".  5  temp). 

i f ( strcmp < stemp.  "he  1 p "  )  ==  0  !  !  s trcmp < s temp.  "HEl P" >==0 
{  h  e  1  p  ( 1  /  <  ; 
goto  qi7; 

> 

sscanf < stemp,  "%f",  f  lofrac+f  lao  >; 

while  (  f  1  Gf  rac  C  f  lag  1  <00  :  I  f 3 ofrac [ f 1 aj 3  >  1  <•) 

<  pr i n t f ( ”\nP 1 ea s e  keep  it  between  0  and  1  "  > ; 
printff"  Try  again 
scanf  i  ”7.f",  f 1  of rac+f 3a„ > ; 

> 

> 

r+f lag. 


/*  go  through  equation-adding  routine  for  each  node  with  dividej  outflow  */ 


1  f  >  f  3  ag  1 ) 

{ 

for  <i  =  l;  i<=flag-l;  i++, 

< 

-r+detrmsz; 

/*  temp  is  the  variable  number  correspondina  to 
this  flow  path  (from  noden  to  tonodeCl)  */ 

temp  =  matr  i  x  [  <  tonodeC  i  3  >  3  Cnodenl, 


determCdetrmsi ICtemp 3  =  (10  -<  f  3  ofrac  [  1  3  >  >; 
for<ii  =  l:  i  i<=<  ary  s  z  t  +  1  >;  ii++> 

<  if  <  lira  tr  i  x  C  i  i  3  Cnodenl  !«  0  )  W  I  ii  '=  tonodeti"i)  > 

determC d e trms z 3 [ < ma t 1 i x C i i 3 [nod en3 ) 3  *  (-1  0  *  ( f 1  of r ac C 1 3 ) ) . 

> 

> 


> 

*  print  out  determinant  for  debug  purposes 
for  <i*l.  1  <=  detrmszi  i++> 

<  printf<"\n")i 

for  (j=l;  j<=detrmsz;  j-r) 
f  printf  <  "7,-7  3f".  d  etermC  i  3  C  j  3  ) ; 

>, 

printf  (  "7.10.  3f",rhsCi3>; 

>. 

/* 


*/ 
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send  determinant  to  routine  solvsys  for  solution 


*/ 

so l * sy  s  <  ) . 

/  *  )}>)>))>>>))>)>)>) i )>)))))))(<'<<(((<<<<(((<((<<''<<<<<  < 

print  out  the  results! 

do  this  both  to  the  standard  \ terminal )  output  a  z 
to  a  file  called  RESULTS 

333333333333333333J333333333CC,::CCCCCCCCCCCCCCCC":CCCCi:C  */ 

/*  open  or  create  the  results  file  */ 

fd  *  op  € n ( "RESULTS” »  1); 

if  (fd  <  0  )  fd  *  creati "RESULTS"*  0644 ) , 

1  see  k  <  f  d .  OL,  2 ) ;  /*  write  the  new  output  at  the  end  of  zt\*  file 

If  an  overwrite  is  preferred/  change  this 
statement  to  1  seek ( f d /  OL, 0 > #  " / 

str c py < 1 inebuf #  "  NEXT  RUN. "); 

linebufCOl  =  *\n'»  linebufC13  =  '\n',  linebufC113  =  \n‘; 

o>rite(fd  /  Imebuf  /  12)# 

for(i=0;  i<»79»  1++)  linebufCi3  *  '\0'; 

spr  intf  (  1  inebuf  #  "\nConstituent:  V,s  Production  level  %2.  Of  ucnstt#  o  1  )  > 
wnte(fd»  linebuf/81)# 

/•*  print  out  pipeflows  between  nodes  */ 

nvar  =  0. 

for  (row=l.  row  O  argsze#  row++) 

<  for  <col-l#  col  O  arysze#  col++) 

<  if  (matr  1 x  Crow!  Ccol  3  !  =  0  ) 

{nvar  =  nvar  r  1# 

pr  i ntf  (  ” \n\n%s  from  V.s  to  V.s  =  7.7  2f"/ 

ucnstt#  nodeCcolI.  ndname*  nodeCrow3.  ndname.  vcrCfnetmCrs^hcoi  j])- 

for  (i=0»  i  079;  i++)  linebufCi3  =  '\0  / 

spr  intf  (  1  inebuf  >  "\n\n7.s  from  */.s  to  V.s  =  7.7  2f"- 

ucnstt#  nodeCcoll  ndname*  nodeCrow]  ndname. varCmatri x i r o*  3  i c  o 1 1 3 / * 
write(fd#  linebuf/81); 


/*  Save  the  water  inflows  to  nodes  */ 
if  ( strcmp ( ucnstt/ "water”  )  ==  0) 

<  nodeCrow].  watrin  *  nodeCrow3.  watrin  +  var Cmatr i x C-ow3C c ol 1 3 . 
f lowsCmatri xCrow3tco) 33  =  var Cmatr i x Cr ow3 C c o l 3 3 . 


Calculate  and  print  out  the  concentration  of  the  c on* t i r jfn t 
The  array  “flows"  stores  the  water  flow  rate  in  ea< h  pipe 
The  "else"  routine  below  calculates  concentration  as  fl**# 
(cubic  meters/time)  divided  by  pollutant  flow  (kg/cime)  times 


a  conversion  factor  of  IOvju  to  make  mg/1. 


else 

< 

if  < f lowsCmatn xirowltcol 33==0)  tone  =  0; 

else  cone  =  <  var  tmatr  i  x  Cr  ouil  [  c  ol  3  3  /  f  1  ou»s  C  rr.a  tr  i  >  r  r  ouil  C  c  o  1 1 3  ) 
*  1000. 0;  /*  convert  kg/cubic  meter  to  rg/1  */ 

pr  int f  <  "NnConc  en tra t  i  on  =  Xg  mg  / 1  “ ,  c  or. c  ) ; 
for  (i=0;  i  <=79;  i++>  linebufCil  =  '\0'< 

spr  int  f  ( 1  ineb uf  i  "\nConc entra 1 1  on  =  V.g  mg / 1  \n “ ,  c  one  ) . 
wr i t  e  <  f  d ■ linebof. Bl); 

> 


/*  print  out  flows  to  outfall  */ 

for  (col  =  1;  col  <=  argsze;  col++) 

<  i f (matr i x [argsz e+1 3 t ccl 3  !«  0  ) 

{  nvar  =  nvar  +  i, 

pr intf < "\n\nXs  from  Xs  to  outfall  =  X7.  2f”, 
uenstt.  nodeCcoil.  ndname.  varimatrixCargsze+llCcoljj); 

for  (  i=0;  i<=79;  i*+)  linebufCil  =  '\0'. 

spr intf ( 1 inebuf . "\n\nXs  from  Xs  to  outfall  =  X7.  2f  ". 
uenstt. nodeCcoij  ndname. varCmatrixCargsze+13lcol3j); 
wr  i  te  <  f  d.  linebuf.81); 


if  < stremp < uenstt. “water "  )  =«•  C> 

f lowsCma tr i x Carg sze+1 3Cco 1 3 3  =  var  Cma tr  i  x  C ar g s  z  en  *  E c  o  1  3  3 ; 


/  * 

Calculate  and  print  out  the  concentration  of  the  constituent 

*/ 


else 

{ 

if  < f lows tmatr i x Carg s z e+1 3 Cc ol 3 3=«9)  cone  =  0. 
else  cone  *  <  var tmatri x Carg sz e+1 3 t c oi j 3  / 
f 1  out Cmatr i x Cargsz i* 1 3Cc ol  3 3  ) 

*  1000.0;  /*  convert  kg/cubic  meter  to  mg/1  */ 

printf ( "\nConcentration  =  Xg  mg/1". cone); 
for  ( i=0.  i  <=79;  i++>  linebufti3  =  '\0'. 
spr  intf  <  1  inebuf .  11  NnConcentrat  i  on  =  Xg  mg/l".conc) 
ur  i  t  e  (  f  d .  linebuf.Sl); 
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/*  print  inflows  to  system  where  there  is  a  flow  requirement  * 

pr l n  t f  < " \n  "  ) . 

for  <i  =  nvar  +  1,  1  <=  detrmsi;  l-*--*-) 

t 

1  f  (  var  C 1 1  1  =  0.0' 

printf  <  "\nlnf  low  to  network  at  7. s  =  7.7  2f\n 
nodeCvarlblCi33  ndn ame,  var Ci 3 ); 

forfcol  =  1;  col<=79;  :ol++)  linebufCcol3  =  ‘\0’ 
spr  intf  <  1  inebuf .  "\n\nInflow  to  network  at  /.s  *  7.7  2f" 
nodeC  var  1  b  1  C  i  3  3.  ndname,  varCi3), 
write(fd,  linebuf,81)> 


else 

{ 

varCil  =  -1  *  varCil; 

printf ( "\nExcess  flow  at  7.s  «  7.7 .  2f\n“. 
node  C  var  lb  1  Ci  33.  ndname,  var  Ci  3  )  ; 

for(col  =  1;  col<”79;  coi++)  1  inebuf Ccoi j  =»  '\C/*i 
spr  intf  <  1  inebuf ,  "\n\nEkce*s  flow  at  V.%  *  \7.  2f 
nodeC  var  lb  1  C  i  3  3.  ndname,  var  C  i  3  ) ; 
wr i te ( f d , linebuf, B1 ), 


unteifd.  '  \n  ' ,  1  ) , 

c 1 ose i f d  ) , 

Use  the  " inf  in i te 1 y-abusab 1 e  goto”  to  allow  rerouting 
of  flows  and  another  run  of  the  program 


ql8  printf ( "\n\nDo  you  want  another  run?  y  or  n  (h/a)  "); 
scanf  (  "7.  s  ",  query); 

if  < strcmp (query,  "help"  )  »«■  0  !!  strcmp  (  query,  "HELP"  )  mm  0> 

{  he  Ip ( 18) ; 
goto  qlS; 

> 

i*  (queryC03  'V'  !!  aueryCOl  “»  y')  goto  recycle; 

>  /*  end  of  bildsys  */ 
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hinclude  <stdio.  h>; 


extern  struct  a  {char  aa£203;  }; 
extern  struct  (char  ndname£20J. 

struct  a  ndtoCdil; 
double  watrin; 

>  node! 30 1 . 

extern  int  matr  l  x  C303C303.  arysie; 
chgmtx (  )  < 

/*  |  i:  |  i  ;  ;  it;;  | 


function  chgmtx  to  make  changes  in  the  flow  connection 
matrix,  allowing  the  correction  of  errors  or  tryinf 
alternate  recycling  schemes 

The  important  variables  are  external/  as  defined  elsewhere 


char  chgnode£20],  elmnodeC201,  addnodeC201; 
int  i,  chgno. 

start  pr int f ( "\n\nWh i ch  node  do  you  want  to  change7  (h/a> 
scanf  (  "V.s" .  chgnode  > . 
i f ' c hgnod eCOl  ==  '0')  return; 

i f < s trcmp ( chgnode, “help " )  ■■  0  !!  strcmp ( chgnode. "HEL P" >  =-  0) 

< 

h  e  1  p  (  7  ) , 

goto  start. 

g. 

/*  find  the  number  of  the  node  to  he  changed  */ 

for  (1  =  1.  l  <=  arysie,  !•*■♦) 

{  i f ( strcmp < chgnode. nodeC i 1  ndnam« )  ==  0  ) 

{  chgno  =  i. 
break; 

> 

> 

if  (l  i  arysie)  <pr int f < "\nNo  Such  node  as  X* " . c h gnoo e ) . 
goto  start; 

> 


print  out  where  flow  ones  from  chgnode 


eight  for  (i  =  l,  i  -:=  arysie;  1++) 

<  i f < ma tr i x C i J t c hgno 1  '  =  0) 

pr int f ( " XnF 1 ow  goes  from  X*  to  X\" , c hgnode. nodeC l 1  ndname). 

> 

l f ( ma tr i x C ar y s i e+ 1 1 C c h gno 1  '*  0) 
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pr- 1  n  t  f  i  " \nFl  ou;  goes  from  V.s  to  outf  ai  1 " ,  c  hgnode  ;  ; 

uih  1 1  e  (  e  1  mnod e C03  !=  'O') 

< 

pr  i  n  t  f  (  " \n\nWh  i  c  h  node  do  unu  uiisft  to  eliminate  flow  to'  (h/a>  ") 
sc anf '  "V.s",  elmnode ) ; 

l  f  <  s  tr  cmp  (  e  lmnode.  "h  el  p  "  >  “*•  0  !!  s  tr  cmp  (  e  lmnod  e.  "HELP  "  )  »=  0) 

{ 

help  (8), 
goto  eight; 


l f ( s tr cmp ( e Imnode. "outf al 1 “ )  ==  0)  <  matr i » Car gs ze+1 j Cc hgno j  -  0. 

i  =  0;  /*  break  out  of  loop 

pnntf  ( "\ndDne" ); 

> 

else 

{ 

for  ii=l;  i  <=  arysze;  l+r) 

{  if  <  strcmp  <  elmnode.  nodsC  i  3.  ndname )  ==  0) 

{  matr i » C l 3 Cchgno j  =  0; 
printf (“\ndone"); 
break ; 

> 

} 

> 

if  ii  >  arysze  ?<&  elmnodeCO]  !*  'O') 
pr  intf  <  "  XnFloie  does  not  go  to  a*” .  elmnode ) ; 

> 

wh 1 1 e < ad dnod e 103  ’=  '0') 

< 

nine  pr  int  f  (  " \n\nWh  ich  node  do  you  wish  to  add  fiow  to"*  <h/a’  ") 
scanf  ( ‘"/is" .  addnode). 

l  f  <  strcmp  <  addnode.  "he  1  p  "  )  ==  0  !!  s  tr  cmp  (  ad  dnod  e.  "HELP  "  >  ==  '•! 

< 

h e 1 p  <9 ) ; 
goto  nine; 

} 

i f ( strcmp ( addnode. "outf al 1 ” )  ==0 )  <  ma tr i < i ary s 2 t* 1 3 C ch gno 3  -  1. 

1  =  0. 

pr  int  f  (  "Xndarne"  > . 

} 

else 

< 

for  (i=l,  1  <»  arysze.  i+r> 

<  if  < strcmp ( addnode. nodeC i 3  ndname >  =•  0) 

<  matr  i  %  C  i  3  C  c  h  ono  j  =  1; 
pr intf ( "\ndone" ) . 
break ; 

> 

> 
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> 


i -f  <  i  >  sri^sie  !<S<  addnod»C03  '*  0'  ) 

pr  intf  (  "NnNo  such  node  a?  ’/.s'\  addnode  ) , 

} 

elmnodeC03  =  addnobef03  » 

goto  start; 

> 


(•include  <stdio.  h>, 

routine  findmdl  to  look  in  a  (ui  of  existing  process 
models  and  read  it  if  it  exists  Also  allows  writing 
over  existing  models. 

Each  model  in  the  file  has  5  lines:  the  first  is  a  title  iin- 
with  an  x  in  the  first  space,  and  4  lines  of  altc-rnating 
x  and  g  points  t  represents  the  production  level,  from 
0  to  100.  V  represents  the  water  demanded  or  mass  of  waste 
chemical  produced. 

Each  line  must  have  80  char,  (not  including  the  newline) 

1 1 1 1 1 1 11 1 1  U  11 1 11 1  111  1 1 1 1 1 1 1 1 1 11  111  1 1 1 11 1 1 1 1 11 '.!  11 1 U 11 1  «/ 


fmdmdKwlabel.  wmd  1.  prolvl.  chngs) 
char  wlabelC203, wmdlC203; 
float  prdlvl. 
int  chngs; 

< 

char  1 inebuf CB1 3.  numbufC93.  fletest; 
char  testC33.  labelC203.  quergC103; 
int  fd,  i,  j ,  nread,  line,  ihi,  ilow; 
float  x  C203,  g  1 20  3 ; 
double  qn, 


/*  initialize  «  and  g  »/ 

for(i»0;  i  <■  19,  <»Ci3  »  0  0;  gti3  ■  0  0;  > 


fletest  =  'n 


/•  fittest  keeps  track  of  what  kind 

•if  writing  is  to  be  done  at  eno 'of  prog 
'r '  *  replace  old  model 

'n'  •  write  new  mooel  in  existing  fils 
'o'  =  use  existing  model,  r.o  overwrite 
c'  •  create  new  file,  writ*  in  ir.udei 


The  variable  ''chngs"  ij»  equal  to  1  i  ' 
chans  in  the  model  are  to  be  anowed  and 
is  0  if  no  changes  allowed 

*/ 
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to  -  oocn  .  uimd  i .  2  »  . 

it  •  f d  '  0) 

{  pnntf ("\nNo  filt  ot  nodtls  found,  file  creat«d\n" / . 
fd  =  c  rea  t  (uimd  1 , 0644 ), 
f  latest  =  'c  '  ■ 

> 

*  If  file  of  models  eiists.  print  out  the  names  of  the  models  * > 
else  i  lsee«;fd.OL.O). 


uihilc((nr«ad  *  r**d<fd»ain«buf/8l)) 

'«  0) 

\ 

if  (linebufCOJ  **  ' x ') 

* 

for  <  i*=0. 

l  <=  IV,  i+r) 

labtl t i 3 

i f < strcmp ( label * wlabel )  «  0) 

*  1 inebuf C  i*l  J 

fietest  «  o', 

/*  desired  model  is  note  found 
reed  it  */ 

ftr  (line  «  0,  line  <»  3,  lin**) 

{ 

read  <  fd.  1  ir.ebuf  ■  81 ) . 
for  ( i  =  0.  i  <*4.  1++) 

{  for(j  =*  0;  j  <=  7,  j-*-+) 

{numbufljl  «  1  inebuf C ( 16*i )  jj, 

> 

numbufin]  «  'SO'. 

sscenf  (numbuf .  "%8f  ",  »+-<l  ’.ne*3  +  i )  )  i 
for<j«0;  j  <*  7,  j++) 

fnumbufCjl  *  1  inebuf t ( 16*i >  *  j  *  o3, 

> 

numbuf  ml  =  'SO'. 

sscenf  <  numbuf ,  "7.Sf".  g*(  1  ine*5  +  ill. 

> 

> 

break, 

> 

> 

> 

/*  test  if  requested  model  not  in  file  */ 
if  If 1 etest  1  *  'o') 

<  pnntf("\nNo  model  celled  7.20s  in  f  1 1  e" .  ui  ab  e  1 ) ; 
goto  build; 

> 

*  if  r.o  changes  are  to  be  made  go  to  end  of  program  */ 
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if  i  chngs  ==  0  $<&  fletest  =  «  'o')  goto  quit. 


/*  plot  out  the  model,  if  ucer  wants  to  change  it. 

or  if  no  model  was  on  the  file,  build  a  new  one  */ 

graph  <  x .  y , 20. wmd 1 ) ; 

pr intf ( " \n  PRODUCTION  LEVEL" ). 

printf("\nDo  you  want  to  change  this  model"  y  or  n  ”), 
sc  anf  ( ’"/.s".  test ) ; 

if(testC01  '=  '  Y '  !<!■<  t  e  s  t  [  0 1  '=  '  y  ' )  goto  a  u  1 1 ; 
fletest  =  '  r  ' ; 


/*  buiia  new  model  */ 

build  testCO]  =  'y 
while  (testCOl  ==  'y') 

< 

for  (  i=0;  i  <=19,  i++)  <xCil  =  0  0;  ytil  =0  0;  labeltll  ”  '  > 

princf("\n\n  3UILD  A  NEW  MODEL  for  V.s\n",  wlabel  ); 

strcpyi'  label,  wlabel); 

1 f ( strcmp < wmd 1 , "water " )  ==  0) 

pr i nt f < " \n Input  water  used  as  a  function  of  production  levelin"). 
else  printfl  "\nlnput  ‘/.s  generated  as  a  function  of  production  levelVn", 
wmd 1 ) ; 

p  r  int  f  <  " \nFor  a  max  of  20  points,  put  in  prod,  level.  ">; 
printf("up  to  100.  \nnot  including  lero,  and  7.s  usa  J/generated " .  wmd  1  ) , 

i f ( strcmp ( wmd 1, “water " )  ==  0; 

pr  intf  ( " \n\nUse  units  of  cubic  meters/day  or  cubic  meters/batch  f9r  waterin''  ), 
else  pr  intf  <  "\n\nUse  units  of  kg/day  or  kg/batch  for  V.s\n".  wmd  1  ) , 
printf>“<for  help  with  units  type  'help')\n"). 


ql5  for <1=0;  i<=19;  i++) 

{  pr  intf  <  "NnProd  leveiCXdl?  (h/a>  ",  i  +  1); 
sc  anf  <  "  V.s".  query  ) . 

l f < strcmp < queru . "hel p " )  ==  0  { !  strcmp ( query. "HELP" )  ==  O' 

{  help  <  15), 
goto  qlS; 

> 

sscanfiquery,  "*/.f",  x  +  i  ) , 
if(iti)  <«  0.0)  { 

pr intf < ”\nPl ease,  no  negative  or  iero  values"), 
goto  build; 

> 

pr  intf  (  "\n7.s?  ",wmdl>; 
scenfCXf*.  y  +  i)i 

if  (xCil  >=  100)  break;  /*  stop  at  1007.  capacity  •»/ 
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> 

graph  ,  t,  y .  20,  wmd  1  ) , 

printf(”\n  PRODUCTION  LEVEL"); 

printfi"\nDo  you  want  to  change  this  model?  y  or  n  " > ; 
scanf  (  "7.s" ,  test ) . 

if(testC01  ==  'Y')  testCOl  =  'y ', 


quit  if  <fletest  '=  '  c  ' ) 

/*  position  file  for  writing  */ 

/*  if  replacing  old  model,  move  back  up  5  lines<405  char)  */ 
<  if  (fletest  ==  'r  <  lseek(fd.  -405L,  1); 


/*  if  writing  new  model  in  file,  move  to  end  of  file  >/ 
if  (fletest  mm  'n  t  lseek(fd,  OL,  2), 

/*  if  creating  new  file,  remind  *■  J 

if  (fletest  ==  'c')  lseek(fd,  OL,  0); 

/*  write  out  the  model,  title  line  and  then  4  data  lines  */ 

I  inepuf  103  =  '  x  '; 

for  < 1  =0,  1  <»  19;  i++)  linebufti  +  U  -  labelCU; 

for  <1  =  21,  i  <*  79;  i++)  linebuf Cil  «  ' 

I I  neb  uf C80 J  =  '\n'; 

write(fd,  linebuf,  91); 


for  (line  =  0;  line  <=  3;  iine++> 

<  for  <i=0;  i  <=  4;  i++) 

<  sprintf  (numbuf.  ”7.S  le",  *1  ( line*5>  il); 
for  <  j=0;  j  <=  7;  j**  ) 

1  inebuf c i*18  +  jl  »  numbuf!  jl, 
sprintf  (numbuf  "7,8  le",  yC(line*5)  *  1J>, 
for  <j»0;  j  <=  7;  j+t) 

1 inebuf C ( i*16)  +j  =  numbuftjl; 

> 

linebufCSOl  »  '\n'; 
writeifd,  linebuf, 81); 

> 

>  /*  end  of  write  routine  */ 

c 1  os e  <  f d  > . 


/*  *««»#*#*#*###  interpolate  value  of  water  demand  ######*## ---###  »/ 


ilow  =  0,  ihi  =  l, 

if<  prdlvl  •;  »t03)  qn  =  gCd;  /*  Assume  that  if  «  less  than 

min  in  graph,  y  stays  at  ymin 
Don't  assume  water  use  g^es  to 
tero  ac  zero  production  */ 
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e  i  s  e 


whiieiprdivl  >  xCihil  &&  ihi  <  20) 

<  ■*■+1  h  i  ; 

+•*•11  ou/. 

> 

ifiihi  >=  20)  /*  D"r't  30  off  the  upper  end  or  the  graph  */ 

<  pr  a  n  t  f  (  " \n‘/.f  V,  production  level  too  high  at  V.s;  extend  ,npdel". 
prdl vl, u/labe  1 ) ; 
goto  build- 

> 


qn  =  yCilou/1  + 

(yCihil  -  ytilouO)  *  ((prdlvl  -  xCilou/3)  / 
(  x  C 1  h  1 1  -  xlilouil)); 

> 

return i qn  ) , 

> 


# 


Routine  Graph  to  mane  a  graph  or  a  tu/o-dimensionai  array, 
such  as  a  process  model 


!  +/ 


graph(  x,  y,  npts.  ylabl  ) 
float  x  [201  ■  yC203; 

1  n  t  npts. 
char  ylablC201; 

< 

float  ymax,  graphy[531,  temp,  xmax.  q> 
1  nt  1 ,  i  l  ou-  ihi,  j> 


/*  *■**  first  find  the  max  x  and  y  values  ***  *t 

xmax  =  0.  0, 
yma«  =0  0, 

for  (i=0,  1  <=  npts-li  i++) 

{  if  (  yCil  >  ymax)  ymax  =  u C i 3 1 
if  (  e  C  i  3  >  inn )  xmax  =  xtil; 

> 

/*  generate  50  points  for  the  graph  by  interpolati on  */ 

1 1  oui  =  0.  ihi  »  1; 
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for  !q*  xmax/50  0;  q  <=  inti;  q  »  q+  (xmax/50.  0)  ) 

{ 

j  =  q  *  ( 50.  0  /  xmax ) ; 

iff  q  <  x  C03 )  graphqC  J  3  =  yCOJ;  /*  Assume  that  if  x  la.s  thar, 

min.  in  graph*  u  stays  at  y.r.i  n 
Don't  assume  water  use  goes  to 
zero  at'zero  production  */ 

else 

< 

u/h  i  1  e  (  q>  x  C  i  hi  3  ) 

<  ++ihii 
++i  low; 

> 

graphyCjl  *  yCilowl  + 

! y [ ih  l  3  -  ytilowJ)  *  <(q  -  x(ilow3>  / 

<  xCihi3  -  x  C i low3 ) ); 

> 

> 

/*  noui  print  the  graph  */ 

lhi  =  s tr len < y lab  1 ) ; 

for  ( i=20;  i  >  0;  i  — ) 

{ 

if(i  <=  ihi)  printf  ("\n7.c".  y  lab  1 C ih i-i 3  ) ;  /*  write  tho  y  label  */ 

else  printf<"\n  "); 

if(  i  V.  4  ==  0)  <  temp  =  (ymax  /  20)  *  i; 

printf  ("7.10.  2f  temp); 

> 

else  printf!"  !"); 

for!j  =  l;  j<=50;  j++) 

{  if  (  graphy  C  j  J  <  <  i+1 )  *<  ymax/20.  0)  !<S<  graphyCjl  >=  i  *<umax/20  O’  ) 
printf! "o" );  /*  print  if  it  is  a  point  */  — 

else  printf!"  /*  print  this  if  it  isn't  a  point  *~ 

> 

/*  print  x  axis  and  label  */ 

pr  in t  f  ! 11  \n  !  "  ) ; 

for  !i*l;  i<“5;  i++) 

printf!  ” - !  "  ) ; 

printf! "\n  ); 

for  <i»l;  i<-5;  i++)  {j  -  xm.x  *  (i/5  0); 

printfC"  X3d".  j);  > 

return; 

> 


#  include  Tstdio  h>, 

help(n) 

int  n, 

< 

/*  help  routine  tor  the  ammo  plant  model  */ 

int  test.  l.  oldn, 
char  quer  y  C  3 1 . 

oldn  =  n . 

wh 1 2  e t  1  ) 

< 

pnntfi"\nDo  you  want  general  heip  or"), 

pr:ntf("  help  specific  to  the  question  you  are  on”'  "); 

p  r  i  n  t  f  i  " \nEnt  er  1  or  2  "). 

prmtf  i"\n\n  1=  general\n  2*sp  ec  i  f  i  c \n" ) , 

scanf  (  "V.d  “ ,  S<test ) , 

i f  ( test  ==  1 ) 

{ 

pr  int  f  <  "  \  n\nAnwium  t  i  on  Plant  Process  Hodel\n"); 

printf < "\n  This  is  a  program  to  model  the  flow  of  water  and  water-borne"), 
pnntf (  ”  pollutants  \nThe  following  general  instructions  applu  "1; 
printf (  "\n\nl  All  process  (node)  names  should  contain  20  o,  fewer"), 
printfC  characters  of\nany  type,  hut  without  blanks."); 
pnntfi"  Be  careful  with  spelling!"); 

pr  i  n  tf  <  ” \n\n2.  Yes  or  no  questions  can  de  answered  with  yes-  no-  y-  or  n  " ) . 
pr  int  f  (  " \n\n3.  Whenever  a  prompt  repeats  itself  you  can  tell  it  to  move  on"), 
printf("\nby  typing  '0'  (zero)  .">; 

p r i n t f < " \n\n4  You  must  always  model  water  flows  before  pollutant  flows  ">. 
printf ( "because  treatmentNnmodel s  use  flow  rate  as  a  parameter  Any  time  "), 
printf  ('"you  change  the  pr  oduc  t  i  on \nl  e  vel  or  re— route  flows  you  should  model") 
printff"  water  flows  again. \n\n" ) ; 

printf ("5.  At  any  prompt  where  '(h/a)'  appears,  you  can  answer  'help'  and  get 
pr  int  f  <  "  it  ">. 

pr  in  tf  (  " \n\n!“lore  details  are  available  in  the  user  manual  and  details  arc")- 
printf<“  avai lab  1 e\nhere  for  specific  prompts  ”), 
p r i n t f ( " \n snHi t  '1'  to  c on t i nue\n" ) ; 

sc  an  f  <  "V.  s  "  ,  query  ) , 

pr i n t f ( "Ch oose  one  of  the  following: “); 


pr  in t  f  < 

"\n\n' 

'*); 

p  r  i  n  t  f  ( 

**  \n 

1 

printfC 

”  \n 

2 

printf ( 

"\n 

3 

pr  i n t  f  ( 

”  \n 

4 

pr intf  < 

"  \n 

5 

pr intf ( 

"\n 

6 

pr 1 nt  f  < 

"  \n 

7 

pr  i n tf  C 

"\n 

S 

Do  you  want  input  to  come  from  a  file?"); 

What  is  the  name  of  the  file?"); 

What  is  the  name  of  the  neat  nnie""1); 

What  node  does  flow  go  to?">; 

Enter  node  wmch  goes  to  outfall"); 

Do  you  want  to  change  flow  connections''"); 

Which  node  do  you  want  to  change?"); 

Which  node  do  you  wish  to  eliminate  flow  to?"). 


51 


printfi''\n  9  Which  node  do  you  wish  to  add  flow  to?"), 
pnntf  i''\n  10  What  production  level?”). 

jrintf ("\n  11  What  constituent  do  you  want  to  model*"). 

pnntf  i"\n  12  Do  you  want  to  use  the  existing  models*'"). 

pnntf  ( "\n  13  Is  there  a  water  demand  at  this  node?'1); 

printf  i"\n  15  far  help  with  units.”); 

pnntf ("\n  17  Enter  fraction  of  flow  from  a  to  b">. 

printf <"\n  18  uo  you  want  another  run?”); 

pnntf ("\n  19  is  the  constituent  generated  hare?"). 

printf ( "\n  24  Is  this  a  treatment  process  where  constituent  is"), 

printft"  removed?"), 

pnntf  ("in  98  for  intr  epretation  of  results"), 
printf i "\n  99  to  quit  HELP  and  return  to  program"). 

pr .nt f i " \n\nWh i c h  one  do  you  want?  "); 

pr  intf  <  "  ( you  entered  'help'  from  no.  Xd  >  ",  oldn., 

scanf  i  "V.d  ",  !<n )  ; 


if  vn==  99)  return; 

else  if  < n==l ) 

<  printf < "\nHelp  for  question  1:  \n"), 

printf ("Do  you  want  process  configuration  input  to  come  from  a  file?"), 
p r intf ( "\n\nl f  you  have  a  file  with  all  the  node  nam^s,  where  flow"); 
pr  int  f  ( "\ngoes  to  for  each  one,  and  the  names  of  nodes  WFTich  drain\n"), 
printf("to  the  out  fall,  yiu  can  use  it  instead  of  answer  >nj\n':  > ; 
p r intf ( "quest ions  3  through  5.  This  is  very  useful  when  running  a"), 
printf ( "Nnmodel  several  times.  If  you  answer  yes,  the  next  question"), 
printf < "\nwi 1 1  ask  you  the  name  of  the  input  file  "), 

> 

else  if  <  n==2 ) 

{  pr i n t f ( "\nhel p  for  question  2  What  is  the  name  of  the  file?"), 

p r i n t f < " \n\n I f  you  told  the  program  that  the  process  configuration"); 
p  r  m  t  f  < " \nwas  in  a  File,  uc.u  now  have  to  tell  it  cne  naipr  of  the"), 
printf ("  file  \nlf  the  file  is  not  in  uour  directory-  be  sure  to"), 
pr  int  f  (  "\ninc  lude  the  path  to  the  file  ’’), 

> 

else  if  (n==3) 

{  printf("\n  help  for  question  3.  What  is  the  name  of  the  next  node?"); 
pr i nt f ( "\n\nYou  now  have  to  enter  the  names  of  all  the  nodes"), 
pr int f <"( proc esses ) \nin  your  model,  one  at  a  time  Each  name  can"); 

pr i n t f ( " Vnhave  up  to  19  characters,  with  no  blanks  The  program"), 

printf("\nis  presently  dimensioned  for  up  to  30  nodes 
pr i nt f <“ \nAf ter  all  nodes  are  entered,  answer  this  quescTnn  with"); 
printf("\na  'O'  (zero)."); 

> 

else  if  ( n««4 ) 

{  printf ( "NnHelp  for  question  4  What  node  does  flow  go  t«?">; 

pr i n t f ( " \n\nFor  the  node  that  you  just  entered,  enter  the  names"). 
prmtf("\nof  the  nodes  which  flow  goes  to  tou  can  enter  from  Q  ), 
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printf!"  to  5\nnodename*.  b*  tun  to  spell  them  correctly  When"), 
pr i n t f ( " \nth er e  are  no  more  to  enter/  enter  s  O'  (zero'-  After  you" 
printf ( "\nhave  entered  all  the  nodes  and  where  flow  goes  for  e»ch"). 
pnntf("\n  one/  the  program  will  tell  you  where  flow  goes  from  each") 
pr int f ( "one\nso  you  can  check  it  “  ) . 

printf("\nlf  it  isn't  right,  you  will  be  allowed  to  change  flow")/ 
p r i n tf < "\nc onnec t i ons  at  a  later  prompt.  ">; 

pr in tf < "NnDon  ' t  worry  about  flows  to  the  outfall  (or  otherwise  out’’). 
printf(”\nof  the  the  siistdm);  you  will  be  allowed  to  tell  the); 
printf ( "\nprogram  which  nodes  go  to  the  outfall  later  ” 


else  if  (n*=5) 

{  printf ( "\nHelp  for  question  5:  Enter  node  which  goes  to  sink”), 
printf ( "\n\nType  in  the  name  of  any  node  (process)  which  drains  to" 
printf ( “\nthe  outfall  or  otherwise  leaves  the  system  'This  premet" 
printf ( “\nrepeats.  so  you  can  have  multiple  outfalls  Enter  '0”). 
printf("  <zero)\nwhen  there  are  no  more. 

> 

else  if  (n==e) 

•£  printf  ( "\nHelp  for  question  6:  Do  you  want  to  cnange  flow  "); 
printf( "connections?"); 

printf ( "\n\nlf  you  want  to  change  where  flew  goes  to  from  any  of  ! . 
printf ( "Snthe  nodes,  answer  yes.  Vou  will  enter  a  routine  which”); 
printf  (  "Nnallows  you  to  do  this."); 

> 

else  if  <n"»7) 

< 

printf  ( "\nHelp  for  question  7:  Which  node  do  you  wane  to  change"”’) 
pr  intf  <  "\n\nYou  are  now  in  the  routine  to  change  flow  foutino  ’’ ) ; 
pr intf ( "\nEnter  the  name  of  the  node  (process)  whose  outflow  you"), 
printf  <  "\nwish  to  redirect.”); 

> 

else  if  (n»*8) 

< 

pr  intf  <  “\nHel  p  for  question  8.  Which  node  oo  you  wish  to  eliminate’’ 
printf!"  flow  to?"); 

printf  ( "\n\nYou  have  just  been  told  which  nodes.  If  aiwj.  flow  goes" 
printf! "\nto  from  the  one  you  want  to  change.  If  you  want  co  stop" 
pr intf < "inflow  from  going  to  any  of  these,  type  in  the  name  of*), 
printf ( "\nthat  node.  Enter  '0'  (zero)  when  you  don’t  want  to’’>. 
printf < "\nchange  any  mors  "); 

printf("\nTo  climate  flow  to  the  outfall  er.ter  the  wori  'outfall 

> 

else  if  <n»«9) 

< 

printf ( "\nHelp  for  question  9:  Which  node  do  you  wish  to  add”). 
printf<"  flow  to?\n\n"); 

printf("If  you  want  to  add  a  flow  path  from  this  node  to  some"), 
printf ( "\nother,  enter  the  name  of  that  node  To  add  flow  to"). 
printf<"\n  tne  outfall,  enter  the  word  'outfall'  Enter  ’O'"), 
pr intf < "\n ( zero )  when  you  don't  want  to  add  any  more  )■ 


53 


> 


else  if  (n==10) 

v  pr 1 nt f < " \nHe 1 p  for  question  10  What  production  level'  i , 

printf i "\n\nEnter  the  production  level  you  want  to  model  This')/ 
pr int f < "\nsnoul d  be  expressed  as  a  percent  of  full  production,  ">. 
pnntf("\nso  should  be  between  0  and  100."), 


else  if  ( n==l 1 ) 

{  p r i n t f ( " \nHe 1 p  for  question  11.  What  constituent  do  you  want"); 
printfl”  to  model?\n\n” ). 

printf(“The  program  will  print  out  a  list  of  the  oollutart"), 
pr int f ( " \nc ons t i tuent s  for  which  models  exist,  anoyou  should"), 
p  r  l  nt  f  (  " \r,c  hoose  one  of  them.  For  the  first  run,  and  after  any"), 
printf ("Vncnang#  in  flow  routing  or  production  level,  you  must  )» 
pr i nt f < "\nmode 1  water  flows  first,  because  the  treatment  models"), 
printf < "Snneed  to  know  the  water  flow  rates. ">; 


else  if  <n»»  12) 

{  pr intf ( " NnHe 1 p  for  Question  12  Do  you  want  to  use  existing  "). 
pr i ntf ( "models?" ) ; 

printf ( "\n\nlf  you  answer  yes,  you  will  net  be  allowed  to  look"); 
printf("\nat  and  modify  any  of  the  process  models.  If  this  ip"); 
printf ( "\nnot  your  first  time  through  the  program  and  uou  are  sure") 
printf  ( "\nthat  you  are  satisfied  with  the  models.  this_<i»ill  save"), 
pr  intf  ( "\nsome  time.  Otherwise,  answer  no  and  look  at.  modify,"); 
pr intf < "Nnand  build  models  as  needed  ">, 


else  if  in==13) 

{  pr intf ( "\nHelp  for  ouestion  13  Is  there  a  water  demand'  ). 

printf ( "\n\nt f  this  node  is  one  where  water  is  used  and  where"). 
printf(”\na  water  demand  vs  production  level  mode)  exists  (or"), 
pr intf ( "Vnshould  exist),  answer  yes  "). 


else  if  ( n== 1 5 ! 

< 

printf < ”\nHelp  for  units  \n\n"), 

printf("To  get  correct  c oncentrat l on  calculations,  all  flows  must"', 

pr  intf  (  ”\nbe  in  cubic  meters/time  and  all  pollutant  gg  nerat  io,,\n"  )  ■ 

pr i nt f < "rat es  must  be  in  k i logrem/time  The  time  value  may  be  any"), 

printf  ( "\nconsi  stent  unit,  such  as  day.  or  by  baten  for  batch’.r."' 

pr i ntf ( "processes  The  concentration  will  be  averaged  over  the"), 

pr i ntf < ”\nt ime  period  chosen. \n\nSome  conversion  factors  \n\n 

pr i nt f ( "Cub i t  feet  x  0283  equal  *  cubic  meters  "), 

printf ( "\nLi ters  x  001  equals  cubic  meters  ">; 

printf ( "NnCal Ions  x  00379  equals  cubic  meters  ">. 

printf ("\nPounds  x  .454  equals  kilograms”). 

pr  i  ntf  (  "  NnCub  i  c  feet/sec  x  2447  equals  cubic  metersxde'^ 

pr int f ( "VnMi 1 1 i on  gal/day  x  3786  equals  cupic  meters/d®y  “). 

printf ( "Nncubic  1 ightyeat s/mi 11 i sec ond  x  1  02ex-40  equals  m3/day”  ■■ 
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else  i#  (n»«17) 

<  pr intf ( "NnHelp  for  auestion  17,  Enter  fraction  of  flow,  >. 

pr int f < 11  \nYou  must  tell  the  program  what  fraction  of  the  outflow"), 
pr 1 nt f ( "Nnf r om  the  first  node  goes  to  the  second.  E  n  c  4  r  a  numoer " ' , 
printf < "\nbetween  jero  and  one"), 


else  if  (n==18) 

{  printf ( "\nHelp  for  question  18.  Do  you  want  another  run"), 
printf<"\nlf  you  want  to  run  the  program  again,  type  'y '  ” ) i 
printf < "NnYou  will  be  allowed  to  change  flow  routing,  model"), 
printf ( "\nanother  constituent,  change  the  production  Jovel."); 
pr  intf  ( "\nand  all  kinds  of  fun  stuff.  "); 


else  if  <n==19) 

<  printf C "NnHelp  for  question  19.  Is  constituent  generated  here"); 
printfCNnlf  this  is  a  process  at  which  the  pollutant  enters  the"), 
printf ( "Nnwastewater  stream,  and  at  which  there  is  a  pollutant"), 
printf  ( "\ngeneration  vs.  production  level  model,  answer  yes”); 
pr intf ( "NnAnswer ing  with  a  tells  the  program  that  there  is"), 

printf  ( “Nnneither  a  pollutant  generation  nor  a  treatment  model”  >; 
printf  (  "\nfor  this  node.  “>; 


else  if  <n»*24) 

<  printf < "NnHelp  for  question  24,  Is  this  a  process  where"), 
printf ("\n  the  the  constituent  is  removed"); 

printf < "Nnlf  this  is  a  process  where  the  oollutant  is  removed"); 
pr intf < "\nfrom  the  wastewater  stream,  sucn'as  a  treatment  process"), 
printf < ”\nand  where  treatment  efficiency  vs.  flow  rate'is  to  te"). 
printf < "Nnmodeled.  answer  yes. ”>; 


else  if  <n*=98) 

< 

printf ( "\nHelp  for  interpretation  of  results  \n\n">; 
printf < "There  are  four  different  kinds  of  output  lines: 

printf ( "\n\n '<constituent>  from  <node>  to  <node>  «  _ '"). 

printf ( "\n\nTh i s  line  tells  you  the  flow  of  the  constituent"), 
printf ( "\n< or  water)  in  the  pipe  between  the  two  nodes  Flow"), 
printf < “Nnrates  are  the  same  as  used  in  the  water  use/pollutant"), 
printf  (  "\ngeneration  graphs.  ”). 
printf < "NnNn 'Concentration  «  _ ’\n"i; 

printf < "NnThis  line  tells  you  the  concentration  of  the  constituent"). 
printf<"\nin  the  pipe  when  modeling  pollutant  flows  ") 
printf  <  "\n\n  '  Inf  low  to  network  at  <node>  =  _ 

pr intf < "\nTh i s  line  appears  when  modeling  water  flows  and  tells  ">, 
pr intM "Nnyou  the  rate  at  which  water  enters  the  system  from"); 
printf  ( "Nnwhatever  the  source  is  <well,  intake  from  rtvar,  “), 
printf  <  "\ntreatment  plant,  etc.  ).  ">; 
pr intf <“ \n\n '£» cess  flow  at  <node>  m 

printf ( “NnThis  line  appears  when  there  is  more  water  flowing  iito"), 
printf ( "Nnthe  node  from  upstream  nodes  than  is  needed,  as”), 
printf  <  "Nndetermined  by  the  water  demand  graph 

pr intf ( “NnTh i s  access  water  leaves  the  system  at  this  point  and"). 
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pr  i  n  t  f  (  "\nc  ou  1  d  be  considered  a*  flow  to  tl,e  outfall  . 

pr i ntf ( "\nThe  presence  of  excess  flow  means  that  you  as e  routing 

pr mtf ( "\ntoo  much  water  into  the  node  ">; 

> 


else 

{  pr mtf ( "\nSorry,  you  picked  an  option  which  isn't  available 
} 

pr  in  t  f  (  "\n\nMor  e  help?  y  or  n  ’*). 
scanf  (  "’/.s  ",  query  ) , 

if  (queryCOl  !■  'Y'  &S<  queruCOl  ' “  'y'>  return, 

> 

> 


^include  -.'stdio.  h>; 
extern  int  matrixC30Ku03; 


*  miiiiii>i>iiiiiiiuiiiiiim<''n)n<>nimxiiiuiixiin) 

subroutine  matrxin  to  build  the  connection  matrix  Read  in 
node  labels,  nodes  being  plant  processes  or  treatment 
processes.  Nodes  representing  sources  to  the  system  (withdrawals 
from  the  river  etc.  )  are  added  automat i cal  1 y . 

variables: 

ndname-  an  array  of  the  node  names.  The  node  number  is  assigned 
as  the  position  in  the  array. 

ndto-  an  array  of  the  node  numbers  vs.  the  node  numbers  flew  goes 
to 


mtxsze-  the  number  of  columns  in  the  matrix.  There  is  one  more 
row  than  columns  to  accomodate  the  sink. 


Modified  to  allow  input  from  a  file. 


struct  a  <  char  aaC202;  >i 

struct  <  char  ndname[203( 
struct  a  ndtoC61; 
double  watrini 
>  nodeC30J< 
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matr x in (  ) 


< 


int  2.  x.  iy .  j.  ndno.  mtx  sz  e; 
char  tosinkC203>  testC43; 
char  filenmC203; 

FILE  *fp<  *f open < ) ; 

on s  printf("\nDo  you  want  process  configuration  input  "); 
printf("to  come  from  a  file?  <h/a>  “)i 
scanf  (  "Xs".  test )  ; 

if  <testC03  ==  'H'  ::  test  C03  =-  “■' )  <  help(l); 

goto  one; 

> 

if  (testco:  !■  'Y'  ScS.  testCO]  !*  'y') 

{ 

pr intf ( “ \nProc ess  conf iouration  entry  routine:  put  in  the  name"): 
pr  intf  (  " \nand  inhere  flow  goes  to  for  each  node.  \n\n">; 


for  <j  =  l,  z  <=  30;  !•«•*-) 

{  three:  pr intf ( "\nWhat  is  the  name  of  the  next  node?”); 
printf("  Hit  zero  for  no  more  <h/a)  "); 
scanf  ( '"/.a",  nodrCz 3.  ndnamei; 
pr  intf  (  "Xs".  nodeCi],  ndname); 
if  (nodeCi  ].  ndnameCO]  •”  '0')  break; 
if  < strcmp <nodeC z ].  ndname. "help" )  0  ! ! 

strcmp  <  nodeC  I  ].  ndname,  "HELP"  )  ■»  0  ) 

{  help (3) ; 
goto  three; 

> 


> 


< 


> 


for  (i  ■  1;  i  <»  5;  i++> 

four:  pr intf (  'NnWhat  node  does  flout  go  to  from  X  a?  (h/a 

nodeC  z  ].  ndname  ); 

scanf  (  "Xs".  nodeC  z  3.  ndtoCi].  aa); 

pr  int  f  <  "‘/.s  " .  node  t  2  3.  ndtoCi]  aa); 

if  (nodeC  z  3.  nd  toC  i  3.  aaCO]  ■■  '0')  break; 

if  (  strcmp  (nodeC 2  3.  ndtoC  i  3.  aa;  "help"  >  “  0  !1 

strcmp  (nodeC  z  3.  ndtoC  i  3.  aa#  "HELP"  )  ■»  0  ) 

<  help  (4); 
onto  four; 

> 


> 


else  < 

tu to  printf  ( "\nWhat  is  the  name  of  the  file?  (h/a)  "); 
scanf  ("Xs".  filenm); 
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if  (  strcmp  <  f  i  1  enm.  "he!  p  "  )  ==  0  !  .'  s trcmp ( fi 1 enm,  "HELP " )  ==  0  ) 

<  help  < 2) ; 
goto  two; 

> 

while  (<fp  =  f open ( f l lenm,  r" ) )  ==  0) 

{  pr  int  f  <  “ \nF  1 1  e  "is  not  found.  Please  tT  g  ocain  ih/e)  w>filer.m). 
sc  anf  (  "'/.s '' .  filenm); 

if  <  strcmp  (  fi  1  enm,  “help '' >  ==  0  I  !  strcmp i fi lenm, "HELP  >  ==  0  ) 

{  he  1 p (2) i 
goto  two; 

> 

> 


for  (  z  =  1 ;  ^z  <=  30;  z++) 

<  f  scanf'l  f  p.  !'/.s".  nodeC  z  3.  ndname  ) , 

if  (nodeC  z  J.  ndnameCul  --  '0')  break; 

for  (i“l;  i  O  5;  i++) 

<  f  scanf  i  fp,  "Xs",  node-Cz  3.  ndtoC  i  3  aa); 

if  (nodeC  1 3.  ndtoC  i  3.  aaCOl  »=  '0')  break; 

> 

> 

> 

mtxsze  =  (i-l); 


/*  The  node  names  and  where  flow  goes  to  have  now  been  read  in. 
Now  construct  the  matrix.  The  column  number  represents 
the  node  flow  comes  from  and  the  row  number  represents  the 
node  flow  goes  to 


for  (i  =  l;  i  <=mtxsze;  i++) 

{  for  <j  =  li  j  <=  6;  j++) 

{  if  <  nodeC  i  3.  ndtoC  j  3.  aaCQ3  ==  '0'  )  break; 

/*  look  up  the  node  with  a  node  no.  equal  to  ndto 

(find  the  node  no.  flow  goes  to.  give  its  name)  */ 

for  (ig  =  1;  ig  <=  mtxsze;  ig++) 

<  if  (  strcmp  (nodeC  i  3.  ndtoC  j  3.  aa.  nodal  ig  3.  ndname  )  ==  0  ) 

<  matri  i  C  ig  3  C  i  3  »  1;  break; 

> 

> 

if  ( ig  >  mtxsze) 

{ 

printfC’r.o  node  called  Jis".  nodeC  i  3.  ndtoC  j3.  aa); 

pr  intf  <  " \r.  Input  error  occurred-  program  abortedNnhnYou  lpse\n\n">, 
exit(  >, 

> 
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> 

> 

/*  add  flows  to  sink  or  outfall  */ 
if(testC03  '«  'V'  «<&  testtO]  •=  'y'> 

< 

for  (  z  =  l;  z  mtiszci  z++> 

<  five:  printf  <  "\nEnter  node  which  goes  to  sink  (h/a)  ’’ )  i 
scanf< tosink); 
if  (tosinkCOJ  =*=  '0')  break; 

1 f ( str cmp ( tosinki "help “ )  *»  O  !!  strcmp ( tosink , "HELP" )  ■■  0) 
<  help ( 5) ; 
goto  five; 

> 

for  (i»l;  iOmtxsze;  i++) 

<  if  ( strcmp < tosink. nodeC i 3. ndname )  =»*■  0  ) 

{  matri > Cmtxsze+1 3 C i 3  «  1; 
break; 

> 

> 

if  (i  >  mtssze)  printf("\nNo  such  node  as  7.%“ ,  to*  ink  1 ; 

> 

> 

else  <  for  (z»l;  z  mtssze;  z++> 

{  fscanfffp,  ::%s",  tosink); 

if  (tosinkCCl  »“  '0')  break; 

for  <i«i;  i  O  mtisze;  i++) 

<  if  <  strcmp  ( tosink.  node!  i  3.  ndname  )  ■«=  0  ) 

<  matri»Cmt»»ze+13ti3  ■  J; 

break; 

> 

> 

> 

f c lose ( f p  )  • 

> 

/*  print  out  labelled  matrix 

jnntf  ("\n  %-20s".  node C  1  J  ndname); 

for  (1=2,  i  <»  mtssze;  i++)  {  pr intf  <  "%-l  5s" ,  nodeC  i  3.  ndname ) ;  > 
for  <i»l,  i  <•  mtssze;  i++) 

<  pr  intf  <  “VnX-tSs".  nodeC  i  3  ndname); 

for  (j“l;  j  <■  mtssze;  j++) 

<  printf  (“•/.-13d».  matr i * C i 3C g3 > ; 

> 

> 

printf < "\nSink  /outfall  "); 
for  <i»l,  i  <»  mtssze;  i++) 

pr  int  f  <  "14-15d" .  matr  i  s  Cmtssze+1  3C  i  3  >;  */ 
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return (mtisz • ) i 


/*  end  matrxin  */ 


♦♦include  <stdio.  h>; 

extern  double  determC 1003 t 1003; 

extern  int  detrmsz; 

extern  double  rhstl003>  varC1003» 

S  Q  1  V  S  IJ  S  (  )  { 

int  i»  iswap.  ii*  row.  coIj 
double  swapo.  mltplr; 
double  c ount ; 


VVW'V^V* 


make  sure  the  diagonal  of  the  determinant  has  no  zeroes- 
this  is  necessaru  to  make  the  solution  routine  work 
First  go  down  trough  the  routs  substituting  them  to 
put  values  in  the  diagonal.  if  this  does  not  fill  in 
the  entire  diagonal,  continue  by  adding  routs. 


sAAAA.-sAA.AA.* 


*/ 


count  =  0.  0; 

for  (i=l;  i  <=  detrmsz.  i++) 

{  if  (determCiKil  *■  0.01 
<  count  «  1.  Ot 

for  (ii  »  l  +  i;  ii  O  detrmsz.  ii*+) 

<  if  ( dctermC i i 1 C i 1  !-  0.0) 

{  for  (iswap  =  It  i  swap  <=  detrmsz.  iswapr+> 
<  sutapo  »  determCi  3C  iswap  It 

determE i 3 C i swap  3  ■  determE  i  i  3 1  isutap  1, 
determE 1 i 3  £  1  swap  3  a  swapo; 

> 

swapo  a  rhsCili 
rhsE i 3  “  rh  si  i  ;  3; 
rhs£ii3  a  swapot 
cent  »  0.  Ot 

break. 


/*  for  diagonal  elements  not  changed  from  zero  by  substitution, 
go  back  through  and  add  an  equation  with  a  non-zero  element 
in  the  appropriate  column. 

•/ 


if  <  count  '»  0.  0) 

{  for  (i»lt  i  <“  detrmsz;  i-1*-) 

{  if  (  determE  i  3E  i  3  «■  0.0) 

{  for  (ii  *  It  ii  <«detrmsz;  ii ♦♦) 
if  (deter, ntii3Ei3  •»  0  0  ) 
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<  for  (iswap«  1;  iswap  <«detrmsi;  iswap*+> 

determC i 3C iswap 3  ”  determC  1  3  C  i  swap  3  *  determri  l  3  C  i swap  3  • 

rhsCil  =  rh*Ci3  +  rhsCiil; 

break, 

> 

> 

> 

> 


/•  Gaussian  Elimination  Routine  (since  the  iterative  one  did'nt 

c onverge ) 

First  eliminate  variables  beloie  the  diagonal  *> 


for  (col  ■  1.  col  <*  detrmsi;  col++> 

{  for (row  «  col  ♦  1;  row  <«  oetrmsi;  row++) 

<  i  f  (  d e termCr  owl  C  c o  1  3  =«  O.O)  continue! 

i f ( determt col  3 Cc ol 3  •«  0.0) 

<  for(isuiap  ■  li  Tswap  <«  detrmsi;  iswap++) 
d e term! c o 1 3 i i swap 3  m  determC c o 1 3 C i swap 3  + 
determCrowlC i swap  3; 
rhsCcol3  *  rhsCcoll  *  rhsCrowl; 

> 

mltplr  »  determCrowl tcol 3  /  determC c ol 3 C c o 1 3; 
for(i=»l;  i<»  detrmsi;  i++) 

<  determCrow3C i 3  ■  determCrowl C i 3  -  (mltplr  *  determCcol 1 C 1 1 ) . 

> 

rhsCrowl  »  rhsirow3  -  (mltplr  *  rhsCcoll); 

> 

/*  debug  print  routine 

printf ( "NnVnDEBUG  PRINTVn”); 

for (row  *  li  row  <“  detrmsi;  row++) 

< 

print f ( “ 
for  ( 

<  print 
p  r  i  n  t  f  (  " 


/*  Now  solve  for  the  variables-  working  up  the  diagonal 
from  the  bottom  right. 

for(row  ■  detrmsi;  row  >*  1;  row — > 

<  varCrowl  «  0.  Ot 

for(col  *  row+1;  col  <“  detrmsi;  col++) 

{  varCrowl  «  varCrowl  -  <  determCrowKcol  3  *  varicoi)>; 

> 

varCrowl  *  (varCrowj  ♦  rhsCrowl)  /  determCrowlCrowl; 


\  n  “  ) ; 

iswap  =  1;  iswap  <■  detrmsi;  iswap++) 
f  (  ”51-7.  3f  determCrowl!  i  swap  1 );  > 

51-7  3f  rhsCrowl  >; 

#/ 
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print  out  variable  values.  (TinE  SOLUTION') 
for  <i  =  l;  i<=d etrir.s 2 ;  i++) 

pr  1  n  t  f  <  " \nvar  lab  1  e  %d  =  V.7  3f  ",  i,  var£  i  3  ) ; 


*/ 


return.  > 

^include  <  s  t  d 1 o .  h> ; 

routine  findmdl  to  look  in  a  file  of  existing  process 
models  and  read  it  if  it  exists  Also  allows  writing 
over  existing  models 

Each  model  in  the  file  has  5  lines:  the  first  is  a  title  line 
with  an  x  in  the  first  space;  and  4  lines  of  alternating 
x  and  g  points.  x  represents  the  production  level,  from 
0  to  100.  Y  represents  the  water  demanded  or  mass  of  waste 
chemical  produced. 

Each  line  must  have  80  char,  (not  including  the  newline) 

liiuiiiuiiiiiiii iiimimir.riiiniimmiiiiniiiiiu  */ 


tr  trod  1 (wlabel, wmdl, qin, chgs) 
char  wlabelC201.  wmdlCiOl; 
float  qin, 
jnt  chgs. 

\ 

char  iinebufCail.  numbufC9],  labeliiOl.  fletest; 
char  test£3L 
char  querg£91, 

int  fd.  i,  j.  nread.  line.  ihi.  ilow; 
float  xC201.  gC201; 
double  qn; 

/*  initiali2e  x  and  g  */ 

for  ( i=0;  i  <=  19,  i++>  <xCil  =  0.0,  gCil  =  0.  0;  > 


fletest  =  'n'» 


/*  fletest  keeps  tracx  of  what  kind 

•f  writing  is  to  be  done  at  end  o-  prog 
'  =  replace  old  mwdel 
'  =  write  new  mod«l  in  existing  file 
=  use  existing  model,  no  overwrite 
=  create  new  file,  writs  in  model 


The  variable  chgs  is  passed  as  a  1  If 
changes  to  models  are  to  be  allowed;  a  0  if  not 


f  d  =  open  ( wmd  1,2); 
if  (fd  '  0) 
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<  pr  i  n  t  f  (  " Nr.No  file  of  models  found#  file  cr eat« d Nn "  ) , 
fd  =  crea t  ( wmd  1 .  0644 ) ; 
f letest  =  'c  '; 

> 

/*  If  file  of  models  exists,  print  out  the  names  of  the  models  */ 
else  <  l  seek  (  fd.  OL#  0); 


whilettnread  =  read ( fd. I inebuf . 81 ) )  '=  0) 

{ 

if  tlinebufCOl  «=*  'x') 

{ 

fort  i=0;  l  O  i9;  i++> 

labeltil  *  1  inebuf C i *1 3 
if  (  strcmp  ( label,  uilabel )  «  0) 

{ 

f’.etest  «=  'o'; 

/*  desired  model  is  now  found 
read  it  */ 

for  (line  *  0;  line  <«  3.  liner*.' 

{ 

readtfd. 1 inebuf . 81 > ; 
for  ( i  *  0;  i  04;  ir+> 

{  fortj  »  0;  J  O  7;  J*+) 

<numbuf[jl  »  1 inebuf C< I6*i >  *  jJ. 

> 

numbufCBl  «  '\0'; 

sscanf (numbuf.  ”7.8f »  +  <  1  '."e*5  ♦  i)), 
fortj«0;  j  O  7;  j+r) 

{numbufCjl  *  1 inebuf C ( lo*i )  ♦  j  ♦  83; 
} 

numbufCi-O  “  'NO'; 

sscanf (numbuf. ”%Bf"> y+< line*5  +  i>>; 

> 

> 

br  eak; 


> 

> 

> 

/*  test  if  requested  model  not  in  file  */ 
if  ( f letest  !*  'o') 

<  printf  <  "NnNo  model  called  7.20s  in  f  i le".  uilabel ) ; 
goto  build; 

> 

/•  if  no  changes  are  to  be  allowed  go  to  end  of  routine 
if  (chgs  *■  0  S*S<  fletest  •«  'o')  goto  quit; 


/•  plot  out  the  model;  if  user  wants  to  change  it. 


or  if  no  model  ujt  on  the  file,  build  >  new  one  */ 

graph  <  x.  y,  20,  "TREATMENT  EFFICIENCY" > , 
printf("\n  FLOW  RATE"), 

pr 1 n t f < “ inDo  you  uant  to  change  this  model"  y  or  n  ”), 

5  c  a  n  f  (  "  X  s 11 ,  test), 

lf(testCO)  !=  'Y'  &&  testCO)  !=  'y')  goto  quit, 
fletest  =  ' r ' ; 


/*  build  neui  model  */ 

build:  testCOl  =  ' y ' , 
while  (test  CO]  ==  'y') 

< 

for  ( i=0,  l  <=19;  i++)  -CxCi]  =  0  0;  yCi]  =0  0;  > 
printf("inin  BUILD  A  NEW  MODEL  for  Xsin " ,  wlab  e  1  ) , 

for<i=0,  i  <=  19;  i*+)  labelcil  =  '  '; 

/*  pr intf ( "\nPut  in  process  nam*  "); 
scanf  (  "‘/.20s",  label);  */ 

printf ( "input  in  removal  percent  as  a  function  of  flow"); 
pr  intf  (  "\nFor  a  max.  of  20  points,  put  in  flow  rate  ">; 
printf  ("and  X  removal"); 

printf  ( "\n\nUse  units  of  cubic  meters/day  or  cubic  meters/batch  for  flowin'1), 
pnntfC'For  help  with  units  type  'help'in"); 

ql5  for<i=0;  i<=19;  i++) 

{  printf ( "inFlow  rateCXdl?  (h/a)  ",  i  +  l); 
scanf  <  "V.s".  query  ) ; 

i f < stremp ( query, "help " )  ==  O  ! !  str qmp ( query , "HELP" )  ==  o) 

<  help ( 15) ; 
goto  ql5; 

> 

s sc anf  <  query,  "‘/.f  ”,  x  +  i  ) ; 

if  (xCO]  <=  0.0)  <  printf  ( "\nP  lease,  no  zero  or  negative  values'')- 

goto  build; 

} 

if  (xCil  <=  0.0)  break; 


printf  ("\nV.  removal?  ">; 
scanf  (  "Xf  ",  y  +  i); 

> 

graph;  x,  y.20.  "TREATMENT  EFFICIENCY"), 
printf! "in  FLOW  RATE"); 

printf ( "inDo  you  want  to  change  this  model?  y  or  n  "); 
scanf  < ”%s", test ) ; 

if(testCO]  ““  'Y')  testCOI  =  'y'; 


quit  if  (fletest  !■  'o') 

/*  position  file  for  writing  “/ 

/*  if  replacing  old  model,  move  back  up  5  lines(405  char)  •/ 
<  if  (fletest  »=  't  ,  1  see  k  (  f  d .  -405L,  1); 
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/*  if  writing  new  model  in  file,  move  to  end  of  file  */ 
if  ifletest  ==  'n'>  lseekifd.  OL.  2>i 

■*  if  creating  new  file,  remind  »/ 

if  (fletest  ==  'c  )  lseekifd,  OL.  0), 

*  write  out  the  model.  title  line  and  then  4  data  Jines  */ 

1 1  neb  uf  CO  1  = 

for  (  1  ■  0;  i  19;  i*+)  1  lr.ebuf  t  i  +  1  3  -  mlabeiCi3, 

for  ii  =  21.  i  O  79;  i++)  linebufCi3  =  '  '. 

1 1 n  e  b  u  f  C  80  3  *  ' \ n  ' . 

mrite<fd,  linebuf.  81). 


for  (line  =  0;  line  <=  3;  line++> 

<  for  (i=0i  i  <=  4.  i++) 

<  spr  intf  inumbuf .  "7.8.  le" .  *  C  ( 1  ine*5 )  +  i3>; 
for  (  j=0;  j  <=  7;  jff) 
lin#bufCi*16  j3  =  numbuflj3; 
spr  intf  inumbuf .  “7.8.  le" .  yC < line*51  +  i3)> 
for  <  j  =  ji  j  O  7;  jft) 

1 inebuf t ( i*16)  +j  +S3  =  numbufCj3; 

> 

linebuft903  =  An'; 
uintei  fd,  linebuf.  81); 

> 

>  /*  end  of  uirite  routine  */ 

c  1  o  s  e  (  f  d  ) . 


/«►  *##»*##*»«*#«  interpolate  value  of  water  demand  */ 


llorn  =  0.  lhi  »  1, 

if(  qin  <  xC03)  qn  =  yC03; 


/*  Assume  that  if  *  less  than 

min.  in  graph,  g  stay  s  at  unun 
Don't  assume  mater  use  gees  to 
tero  at  tero  production 


else 

i 

uhile(qin  >  «Cihi3  &&  lhi  <  20) 

<  ++ihi; 

++i  low. 

if  < i h i  >«  20 i 

{printf  ( "\nFlom  rat#  of  V.f  too  high  at  Zs.  extend  model, 
qin.  wlabel ); 
goto  build; 

> 

qn  =  gCilowl  ♦ 

(gCihil  -  gCiloml)  *  <<qin  -  xCilowl)  / 

( x  C  i h  i  3  -  «t  i  low3 ) ); 


returniqn) . 

> 
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